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ABSTRACT 

An  often  ignored  aspect  of  electromagnetic  radiation  from 
antennas  is  the  characterization  of  their  near-fields.   A 
computer  program,  Numerical  Electromagnetics  Code  (NEC),  is 
validated  for  accurate  near-field  computations  and  applied  to 
a  model  of  a  broadcast  monopole.   E-  and  H-fields  are  plotted 
as  a  function  of  position  along  the  antenna  for  various 
distances  from  the  surface.   The  fields  are  also  plotted  as 
a  function  of  radial  distance  outward  for  various  heights. 
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I.   INTRODUCTION 

A.   PROBLEM  ENVIRONMENT 

Maxwell's  equations  describe,  for  all  space,  the  Electric 
(E)  and  Magnetic  (H)  Fields  that  are  the  medium  for  trans- 
mitting information  (virtually)  instantaneously  over 
extended  distances.   This  practical  application  of  electro- 
magnetic (EM)  radiation  has  generated  considerable  interest 
in  its  far- field"  characteristics.   The  E  and  K  vectors  can 
be  analytically  determined  when  the  current  (I)  on  the 
radiator  (from  which  the  field  vectors  are  generated)  has 
been  described.   For  the  far-field  region  of  monopoles  and 
dipoles ,  acceptable  results  are  obtained  by  assuming  a 
sinusoidally  varying  current.   However,  when  the  field 
vectors  are  calculated  for  field  points  ar  distances  less 
than  X  /  2tt  (the  near-field  region),  significant  errors  occur 
between  calculated  values  and  empirically  measured  values. 
This  discrepancy  is  not  the  result  of  inaccuracies  in 
Maxwell's  Equations,  but  ean  be  wholly  attributed  to  an 
inaccurate  descrioticn  of  the  current. 


"The  far-field,  for  a  monopole  or  dipole,  is  defined  as  dis- 
tances from  the  radiatior  greater  than  X/2tt,  where  X  is  the 
wavelength.   For  free  space,  \-cif   where  f=the  frequency  at 
which  the  radiator  is  excited,  and  c=the  speed  of  light  in 
free  space. 


In  the  past,  the  inaccuracies  of  the  near  field  calcula- 
tions were  not  a  problem — there  was  no  practical  application 
to  which  the  energy  could  be  applied.   However,  in  recent 
years  there  has  been  an  escalating  interest  in  the  close-in 
problem.   This  interest  has  been  generated  by  an  increasing 
proliferation  of  devices  which  use  Radio  Frequency  (RF)  energy 
in  such  a  manner  rhar  The  user  is  exposed  to  near-field 
radiation  (e.g.  microwave  ovens,  hand-held  walky- talkies , 
wireless  telephones,  RF  sealers,  etc.). 

The  U.S.  Navy  was  one  of  the  first  organizations  to 
express  an  interest  in  these  fields.   The  interest  was  the 
result  of  a  concern  for  the  hazards  these  fields  might  pose 
to  fuels,  ordnance,  other  Ell  operational  and  test  equipment, 
and  to  personnel.   The  hazards  posed  by  fuels,  ordnance, 
and  other  FM  equipment  result  in  overt  effects ;  the  hazards 
posed  to  personnel  are  not  so  obvious. 

With  virtually  unanimous  consensus,  the  literature 
identifies  the  heating  of  body  tissue  as  the  primary  hazard 
to  personnel.   However,  some  distinctly  non-thermal  effects 
have  been  observed.   Some  of  these  effects  [Ref.  1]  include: 

a)  minor  changes  in  blood  properties 

b)  a  "buzz"  heard  by  certain  people  when  exposed  to 
microwave  radiation" 


"The  "buzz"  is  assumed  to  be  a  function  of  the  pulse  repetition 
rate  (?RR)  rather  than  the  carrier  frequency. 


c)  abnormalities  of  the  chromosome  structure 

d)  movement,  orientation,  and  polarization  of  protein 
molecules 

e)  changes  in  the  transport  rate  at  the  blood-brain 
barrier 

f)  comfort  imbalances,  such  as  epigastric  distress, 
emotional  upset,  and  nausea. 

The  significance  of  these  effects  is  nor  yet  understood. 

The  general  public  has  not  been  exposed  to  near-field 
radiation  hazards  of  sufficient  intensity  or  duration  tc  be 
cause  for  alarm;  however,  there  are  certain  vocations  in 
which  the  incumbents  have  an  elevated  probability  of 
exposure  to  this  type  of  field." 

B.   THESIS  STATEMENT,  SCOPE  AMD  LIMITATIONS 

In  order  to  provide  more  accurate  data  for  people 
researching  biological  effects  of  near-field  radiation,  this 
study  will  validate  a  computer  program  which  computes  the 
near  electric,  magnetic,  and  composite  (peak)  fields  of  an 
arbitrary  radiator.   The  program  will  also  be  used  to 
compute  these  near  fields  for  a  model  to  be  described  below. 
The  program,  Numerical  Electromagnetics  Code  (NEC),  was 
developed  at  Lawrence  Livermore  Laboratory  under  the 


"Sailers  at  sea  or  steeplejacks  who  effect  antenna  repairs , 
for  instance . 


sponsorship  of  the  Naval  Ocean  Systems. Center  (MCSC)  and  the 
Air  Force  Weapons  Laboratory.   It  can  be  easily  learned  and 
applied  directly  by  a  researcher  who  has  some  knowledge  of 
electromagnetic  theory. 

The  model  to  be  used  for  The  investigation  will  be  a  75 
meter  high  broadcast  monopole  radiating  above  a  perfectly 
conducting  ground  plane.   The  ground  plane  is  located  in  the 
X-Y  plane;  the  monopole  is  co-axial  with  the  Z-axis.   The 
antenna  will  be  excited  at  its  base  with  a  voltage  sufficient 
to  produce  a  radiated  power  of  100  0  watts.   The  excitation  is 
at  a  frequency  of  1  Megahertz,  which  has  a  wavelength  of  300 
meters;  hence,  the  monopole  is  a  quarter-wavelength  radiator 
operating  at  resonance.   A  diagram  of  the  model  is  shown  in 
Figure  1. 

In  simulating  this  model,  NEC  requires  that  the  antenna 
be  broken  into  virtual  short  straight  segments.   In 
consonance  with  this  requirement,  three  different  configura- 
tions will  be  used  for  computations;  a  model  composed  of  5, 
15,  and  2  5  segments  of  15,  5,  and  3  meters,  respectively. 

Using  cylindrical  coordinates,  the  radial  and  parallel 
electric  fields,  and  the  tangential  magnetic  field  will  be 
computed.   Field  values  will  be  computed  in  .1  meter 
increments  for  points  parallel  to  the  z-axis  at  distances  of 
1*  2,  5,  10,  10  0,  and  10  0  0  radii  (1  radius  =  .3  meters). 


'1  radius  =  surface  of  the  monopole. 


X  *  f  <'  f\ite. 
Ground 

Plane, 
Perfectly 

C  onduc  t  »'ng 


Y 


*1  SHI'S 
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Field  values  will  also  be  calculated  along  the  radial  axis 
(Y-axis)  in  increments  of  .5  meters  for  heights  (h)  of  2h, 
Ih,  l/2h,  l/4h,  and  2/100ths  h.* 

The  near-field  examined  in  this  study  is  the  cylindrical 
volume  bounded  by  a  cylinder  of  one  wavelength  radius 
centered  on  the  Z-axis,  rising  from  the  ground  plane  to  a 
height  twice  that  of  the  monopole. 


This  is  the  approximate  height  of  an  average  mans  upper 
torso . 


11 


TT 


DESCRIPTION  OF  SELECTED  ANALYSES 


A.   CLASSICAL  ANALYSIS 

The  antenna  utilized  in  this  study  is  a  member  of  the 
class  identified  as  'thin-wire'  antennas.   The  criteria 
which  adumbrate  the  thin-wire  approximation  are  [Ref.  2]: 

a)  Transverse  currents  are  negligible  relative  to  axial 
currents 

b)  Circumferential  variations  of  the  axial  current  are 
negligible 

c)  The  current  can  be  represented  by  a  filament  current 
on  the  wire  axis 

d)  Boundary  conditions  on  the  electric  field  need  be 
enforced  only  in  the  axial  direction. 

By  far,  the  largest  subset  of  thin-wire  antennas  is  a 
class  known  as  short  antennas.   Shorn  antennas  are  defined 
as  having  one  dimension  (h)  much  greater  than  zhe    other 
dimensions,  but  h  is  also  much  less  than  a  wavelength, 
typically  h  <  \/8.   Because  of  the  predominance  of  this 
sub-set,  it  is  reasonable  that  it  should  be  used  as  a  model 
for  classical  analysis. 

Based  on  Maxwells  Equations  and  geometric  considerations 
as  shown  in  Figure  2,  the  following  non-zero  complex  scalar 
field  components  can  be  derived: 


>-Y 


Figure    2 


gSl        sine    (li  +    1      _   -lT)e-^r  (1) 

M-tt  r  2  3 

r  3r 


E      =    nlh   cose    (Q    +    1      _      j    )e^Br  (2) 

r         27T  r2         ^7 


K      =    £±   sin9    (li  +    1      +    0)e"jBr  (3) 

$        4tt  r  2 


where:  n  =  the  intrinsic  impedance  of  the  medium,  I  =  antenna 
current,  X  =  wavelength,  S  =  2-rr/X  =  phase  constant,  h  = 
monopole  height,  r  =  distance  of  field  observation  point 
from  monooole  base.   Since  -here  is  no  variation  in  * : 


<f>     9    r 


The  interpretation  of  these  equations  under  classical 
analysis  require  two  additional  assumptions: 

a)  The  distance  to  the  field  observation  point  is  large 
relative  to  the  monopole  height,  and 

b)  An  harmonically  time  varying  current  is  uniform  over 
the  height  of  the  monopole. 

The  first  of  these  assumptions  is  made  to  focus  attention  to 
that  volume  of  space  where  the  lines  of  force  have  detached 
themselves  from  the  antenna  and  are  aporoaching  the  form  of 

14 


a  plane  wave  front.   The  second  assumption  is  to  form  a 
first  approximation  of  the  current  distribution  on  the  wire 
which  is  solvable  in  closed  form.   Extensive  analysis  over 
many  years  has  shown  this  to  be  a  valid  assumption  when  the 
antenna  is  excited  at  a  single  point. 

In  this  study,  we  will  be  interested  in  fields  where  the 
observation  points  are  not  at  distances  considered  large 
relative  to  the  monopole  height.   When  this  is  the  case,  the 
diagram  of  Figure  2  must  be  modified  as  shown  in  Figure  3*-':. 
From  the  geometry  of  this  figure,  the  following  non-zero 
complex  scalar  field  components  can  be  derived: 


jnl     -]6rl     "]3r2  ^BrQ 

E   = r-2.  (S +  ® 2(cos8h) )      (5) 

9     ^T    rx       r2  rQ 


F     J   0  (z-h  e _  +  £lil  ~ 1   2zcos3h  e \ 

*"p     4tt     p     r1       p     r2         o       rQ 


(5) 


H#  =  ^  (e       +  e       -  2CcosBh)e     J)  (7) 


"Equations  (1)  through  (3)  are  expressed  in  spherical  coor- 
dinates while  Equations  (5)  through  (?)  have  been  transformed 
to   cylindrical  coordinates. 
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A  current  of  the  form 


Hz)  =  I^singh  for  0  <_  z  <_  h  (3) 


has  been  inserted  into  Equations  (5)  through  (7). 

Equations  (5)  through  (7)  indicate  that  each  field 
component  is  the  phasor  sum  of  contributions  from  three  point 
sources.   The  first  source  is  located  at  the  top  of  the  mono- 
pole,  a  second  image  source  is  located  in  a  symmetric  position 
when  the  first  source  is  reflected  through  the  ground  plane. 
The  third  source  is  located  at  the  base  of  the  monopole;  its 
relative  amplitude  is  modified  by  the  (-2cos8h)  factor. 

Since  we  have  a  solution  for  computing  the  field  points, 
it  is  now  reasonable  to  determine  the  range  over  which  it  is 
valid.   By  inspection  we  note  that  for  r  =  \/2it  in  Equations 
(1)  through  (3),  the  magnitude  of  each  non-zero  term  in  the 
parenthesis  is  equal.   This  implies  that  r  =  X / 2tt  may  be  a 

transition  point.   Further  investigation  indicates  that  for: 

3  9 

a)  p>X/2tt:  the  1/r   terms  in  E   and  H  ,  and  the  1/r 

term  in  E   are  dominant 
r 

3  9 

b)  r<\/2-n:    the  1/r   terms  in  E,  and  E  ,  and  the  1/r 

9       r 

term  in  H.  are  dominant. 
Since  we  know  that  power  drops  off  as  the  inverse  square  of 
distance,  and  that  the  power  flow  at  a  field  point  can  be 
defined  in  terms  of  its  E-field  as: 


17 


P  =  E2/n  (9) 


2  2 

we  can  extrapolate  from  a)  above  (E  8/n«   k/r  .  where  k  =  a 

constant)  that  the  far-field  begins  at  r  =  X/2n    and  continues 

out  to  infinity,  while  the  near-field  begins  at  the  monopole 

and  continues  out  ~o  r  =  X/2rr.;': 

It  is  interesting  to  note  the   complementary  relationship 

that  is  maintained  by  the  power  in  this  transition  region. 

Instantaneous  power  flow  at  a  point  in  space  is  given  by  the 

complex  Poynting  vector: 

S*  =  1/2  (E  x  H)  (10) 

For  the  short  monopole,  the  complex  Poynting  vector  is 
oriented  in  the  radial  direction  and  given  by: 


? 

— ...   _  n  /  In  w  .  „  v  -  j  S  r *  2  . 8  1    \  *                                        rm 

S"   =  tt   (77—  )(sm9)e  J   )  (— t*-  -  -^-r)r-                (ll) 

r  Sr 


The  real  oart  of  this  equation  gives  the  average  rate  of  flow 

2 
of  Real  power  (<*I/r")  while  the  Imaginary  part  is  proportional 

to  the  difference  between  average  magnetic  energy  density 


"This  argument  is  offered  not  as  a  proof,  but  as  a  reascnabl 
rationale  which  can  be  envoiricallv  verified. 


and  average  electric  energy  density  C«l/r  ).   Interpretation 
of  this  equation  indicates  that  for: 

a)  r>x/2Tr:  Reactive  power  flow  decreases  rapidly  O— r-)  , 

1       r 
Real  power  flow  decreases  not  so  rapidly  («—«-) 

r 

b)  r  =  A/2u:  Average  Real  power  flow  and  average  Reactive 

power  flow  are  equal 

c)  r<X/2ir:  Reactive  power  flow  becomes  much  greater 

5  2 

(«r  )  than  Real  power  flow  (=r  ). 

The  transverse  component  of  the  complex  Poynting  vector: 


S*   =  n(—  e"^6r)2cos9sine(-i(^-T  +  -K) )  9  (12) 

r    8r 


is  completely  Imaginary,  indicating  no  Real  power  flow,  and 

it  falls  off  rapidly  with  increasing  distance  due  to  1/r   and 

5 
1/r   contributions." 

The  deviant  behavior  of  the  fields  encountered  at 

distances  r<X/2-r  is  the  result  of  moving  from  a  region  where 

the  lines  of  force  are  closed  on  themselves  and  approaching  a 

plane  wave  front,  to  a  region  where  the  lines  of  force  still 

maintain  a  strongly  spherical  wave  front  and  some  of  the  lines 

are  still  attached  to  the  antenna.   Even  in  the  simplest  of 


"Also  noted  is  the  fact  that  the  Electric  field  changes  from 
Elliptical  polarization  to  Linear  polarization  in  passing 
from"  r>\/2iT~  to  value  of   r<X/2-ir  when  9  t    0°,  90°,  or  180°. 
For  these  angles  the  polarization  of  the  electric  field  is 
linear  on  both  sides  of  r  =  X / 2tt  . 
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cases,  the  current  distribution  in  conductors  is  complex. 
This  complexity  is  the  result  of  non-uniform  currents  and 
charge  densities,  capacitance  and  inductance  generated  by 
these  non-uniformities ,  and  discontinuities  from  impressed 
forcing  functions . 

It  is  well  established  that  as  the  field  observation  point 
moves  within  the  \/2tt  distance  from  the  antenna,  that  the 
measured  values  deviate  from  those  predicted  by  Equations  (5) 
through  (7).   It  is  also  well  established  that  Maxwells 
equations  are  valid  for  all  space.   Consequently,  the  con- 
clusion is  That  the  sinusoidal  current  distribution  is  not 
sufficiently  accurate  for  determining  field  points  within 
the  near-field  region.   The  problem  then  Is  to  find  a  better 
description  of  the  current  distribution.   This  is  the  domain 
of  contemporary  analysis. 

Before  proceeding,  it  is  necessary  to  reflect  that  the 
analysis  performed  thus  far  has  been  for  a  short  antenna. 
The  model  for  this  study  is  a  resonant  monopole;  when  the 
antenna  is  lengthened  to  resonance  some  interesting 
characteristics  appear.   Primary  among  these  is  the  fact 
that  the  reactive  component  of  the  impedance  vanishes ,  which 
says  the  input  reactance  goes  to  zero.   The  current  on  the 
antenna  and  the  voltage  required  to  place  it  there  are 
determined  solely  by  the  resistance,  the  real  part  of  the 
impedance.   When  a  sinusoidal  current  is  assumed  on  the 
antenna,  the  tangential  field  component  is  given  by: 
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The  (-2cos8h)  term  representing  radiation  from  the  base  of 
the  monopoie  has  disappeared.   The  monopole  appears  as  a 
point  source  radiating  from  its  height,  and  an  image  source 
found  at  the  reflection  through  the  ground  plane. 

Also,  the  resonant  condition  occurs  at  a  length  less 
than  that  expected  from  the  X/M-  relationship.   This  is  the 
result  of  an  "effective"  lengthening  of  the  anienna  by 
capacitance  and  "fringing"  at  the  end-cap  of  the  wire. 
Capacitance  is  caused  by  a  build-up  of  charge  at  the  end  of 
a  wire  of  finite  width.   The  change  in  direction  of  the 
geometry  at  the  end  of  the  wire  also  changes  the  direction 
of  the  lines  of  force  emanating  perpendicular  from  it;  this 
is  the  fringing  effect.   These  factors  cause  resonance  to 
occur  for  an  antenna  length  of  h  =  X/M-  -  5  where  5  is  a 
function  of  the  capacitance  caused  by  the  two  effects 
discussed  above. 

The  short  monopole  development  discussed  above  is  very 
restrictive  and  cannot  be  applied  directly  ~o  resonant 
monopoles.  However,  when  the  resonant  peculiar  characteris- 
tics mentioned  are  taken  into  consideration,  the  resulting 
development  will  follow  the  same  procedure  as  used  for  the 
short  monopole.  For  the  resonant  monopole  development  and 
results,  3ee  Ref.  4. 
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B.   CONTEMPORARY  ANALYSIS 

Schelkunoff  [Ref.  3:  pp.  370-3  74]  identifies  three  methods 
for  obtaining  a  closer  approximation  to  the  actual  current 
distribution  on  an  antenna.   They  are: 

a)  The  Integral  Equation  method 

b)  The  Sweep-Off  method 

c)  The  Mode  Theory  of  Antennas  method. 
Each  are  discussed  below,  in  succession. 

1 .   The  Integral  Equation  Method 

The  assumption  is  made  that  the  curreni  is  distributed 
on  the  surface  of  a  hollow  cylinder  which  is  divided  into 
filaments  of  angular  density  (I(zf)/2ii)  and  angular  width 
(d$').   The  Kernel  Eunction  is  modified  to  allow  for 
variations  in  $  and  is  expressed  as  an  integral  around  the 
circumference  of  the  wire  as : 
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The  field  intensity  for  a  straight  current  filament: 
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is  modified  by  inclusion  of  Equation  (1^),  with  the  boundary 

condition:   (E„  +  E   =0)  enforced  to  give: 
z     z  & 
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This  Integral  Equation  (from  which  the  name  of  The  method 
derives)  is  the  limit  (as  N-*«)  of  a  series  of  equations  of 
the  form: 
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which  is  also  of  the  form 


-EX(z)  =  Z    I(z)  (13) 
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Equation  (15)  is  a  'circuit  equation'  for  an  antenna  and  can 
be  solved,  in  a  manner  analogous  to  the  use  of  Kirchhoffs 
Equations  for  lumped  networks,  to  an  arbitrary  degree  of 
accuracy  through  the  use  of  Equation  (17). 


It  is  instructive  at  this  point  to  pause  and 
consider  how  accurately  our  initial  approximation  of  the  current 
approached  reality.   Integrating  Equation  (15)  by  parts  and 
setting  -E  (z)=0,  as  it  is  on  the  surface  of  the  antenna,  we 
obtain: 


(I(z')v(z,zf)  +  Hz')^-)  -   /  2  (l-L-  +  32I)v(z,zT)  dz'=0 

32    z]      dzT^ 
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This  equation  shows  that  as  the  antenna  radius  approaches 

zero,  the  first  term  of  Equation  (19)  approaches  a  constant 

limit  for  all  z,  except  z  =  z-  or  z  =  z„.   Sufficiently  far 

from  the  ends  this  constant  is  small.   The  v(z,z')  factor  in 

the  second  term  is  infinite  at  z  =  z f  and  large  in  its 
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vicinity;  hence,  the  ( — -*■  +  8  I)  factor  must  be  small. 

dzf- 

Therefore,  as  the  radius  of  the  antenna  approaches  zero, 

2_ 
,dl    2_v         ,  -.,._.    _ ,  ,         , 

( ^-+8  I)  approacnes  zero,  wnicn  implies  I(z)  approaches 

dz'Z 
a  sinusoidal  form. 

2.   The  Sweep-Off  Method 

This  method  starts  with  Equation  (15),  the  tangential 

field  equation  for  a  filament  current,  being  integrated  twice 

by  parts.   Then  an  approximation  to  the  current  distribution 

on  a  transmitting  antenna  of  finite  radius,  assumed  to  be 

given  by: 
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I(z')  =  A  cosBz'  +  3  sinSz'  (20) 

is  substituted  into  the  integrated  form  of  Equation  (15)  to 
yield: 


-j6r2  -JSr1  -jSr2 

E^  =  J [I'(zJ^ -  I'(z.)- +  I(z0)|~(- -) 
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This  equation  is  an  exact  expression  for  the  electric  field 
intensity  parallel  to  a  sinusoidal  current  filament.   The 
antenna  is  then  subjected  to  a  compensating  impressed  field, 
-E  ,  and  the  correction  terms  from  the  first  aDDroximaf ion 
are  computed.   The  name  for  this  method  results  from  the 
concept  of  ' sweeping-of f '  the  surface  of  the  antenna  the 
residual  tangential  electric  intensity  by  applying  an  equal 
and  opposite  intensity. 

It  should  be  noted  here  that  although  Equation  (21)  is  an 
exact  exnression  for  E  ,  that  E   does  not  vanish  along  the 
filament.   This  apparent  dichotomy  results  from  the  erroneous 
assumption  of  uniform  convergence  for  E  .   In  actuality,  the 

convergence  of  E   is  non-uniform  at  o  =  0 ,  where  p  is  the 

°  z 

distance  from  the  filament.   For  an  antenna  radius  r  =  a, 


there  exists  a  remainder  term  I,(z)  in  addition  to  I(z); 

In(z)  is  sufficient  to  cause  E  (a)  =0.   As  a^-0  ,  I_(z)-*-0: 
1  z  1 

conseauentlv,  its  contribution  to  E   for  anv  fixed  value  of 
*■  J  z 

p,  greater  than  a,  will  approach  zero  with  a.   In  the  limit, 
E   is  given  by  Equation  (21)  for  ?>0,  but  at  p  =  0  Equation 
(21)  does  not  hold. 

3 .  Mode  Theory  of  Antennas 

This  method  is  based  on  the  solution  of  Maxwells 
Equations  subject  to  boundary  conditions  at  the  antenna 
surface,  and  the  surface  of  the  source  of  power.   The  method 
consists  of  calculating  modes  of  propagation  consistent  with 
the  boundary  conditions  at  the  lateral  surfaces  of  the 
antenna,  then  the  modes  consistent  with  free-space  propaga- 
tion.  These  modes  are  combined  to  satisfy  all  boundary 
conditions  and  the  current  is  obtained  as  the  sum  of  two 
components:   the  TEM  cr  Principal  wave,  and  the  Complementary 
wave  consisting  of  all  higher  order  modes.   The  principal 
mode  is  governed  by  ordinary  transmission  line  equations 
(in  terms  of  distributed  capacitance  and  inductance). 

4 .  Discussion 

Of  the  methods  discussed  above,  ' Sweeping-Of f '  is  the 
most  elementary  and  consequently,  the  least  precise.   The 
Mode  Theory  of  Antennas  excites  a  current  in  the  conductor 
for  the  TEM  mode,  and  a  current  for  the  complementary  wave, 
which  is  a  sum  of  currents,  one  for  each  of  the  higher  order 
modes  excited.   The  excited  modes  are  strongly  dependent  en 
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the  geometrical  configuration  of  the  system;  hence,  a 
generalized  solution  requires  a  number  of  terms  which  rapidly 
approaches  infinity.   This  overwhelming  problem  of  accounting 
for  all  terms  limits  the  attractiveness  of  this  method. 

The  Integral  Equation  method  appears  most  promising 
for  generating  a  better  approximation  to  the  current.   Many 
different  approaches  to  solving  the  equation  have  been 
proposed  and  virtually  all  of  them  lend  themselves  to 
solution  by  digital  computer.   Three  of  these  approaches  will 
be  discussed  in  the  next  section. 

C.   NUMERICAL  ANALYSIS 

1 .      Hallen's    Approach 

One  of  the  earliest  solutions  to  the  Integral  Equation 
was  the  approach  proposed  by  Hallen,:';  the  basis  of  which  is 
successive  approximations  to  the  current.   Kallen  uses  a 
traditional  Green's  Function  along  with  a  tacit   assumption 
of  a  homogeneously  distributed  current""  over  the  conductor 
surface.   The  'ET  field  inside  the  conductor  (expressed  in 
terms  of  current  and  skin  effect  resistance)  is  equated  to  the 


"Hallen's  approach  is  relatively  long;  the  details  will  not 
be  presented  here.   For  an  outline  of  the  method  and  a  discus- 
sion of  the  important  steps  and  results,  see  [Ref.  S]. 

""Hallen's  tacit  assumption  is  that  the  impressed  voltage  is 
distributed  over  a  finite  length  [Ref.  3]  (which,  for  a  homo- 
geneous conductor,  implies  a  homogeneously  distributed  current) 
If  the  stated  assumption  of  a  driving  potential  discontinuity 
is  strictly  enforced,  Hallen's  method  breaks  down  [Ibid.]. 
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'E'  field  outside  "the  conductor  (expressed  in  terms  of 
vector  potential,  'A').   The  wave  equation  is  solved  as  a  sum 
of  a  complimentary  function  and  a  particular  integral.   A 
constant,  C~,  in  the  solution  is  evaluated  in  terms  of  input 
conditions  at  the  terminals,  and  the  vector  potential.  A,  is 
expressed  in  terms  of  the  antenna  current.   C?  and  A  are 
inserted  in  the  solution  for  the  wave  equation  obtaining  an 
Integral  Equation  in  current,  which  is  Hallen's  Equation.   A 
partial  solution  for  the  current  is  obtained  by  evaluating 
one  of  the  integrals  so  the  current  is  expressed  as  the  sum 
of  several  terms,  some  of  which  also  include  the  current,  I. 
Neglecting  certain  terms  in  I,  an  approximate  (zero  order) 
solution  is  obtained  for  I.   This  value  is  substituted  back 
in  the  current  equation  obtaining  a  first-order  approximation, 
This  process  is  continued  to  yield  successively  higher  orders 
of  solution.   The  process  is  stopped  when  a  sufficient  degree 
of  accuracy  is  obtained,  the  constant  C1  (the  coefficient  for 
the  quadrature  term  of  current)  is  evaluated  and  an 
asymptotic  expansion  for  the  current  is  obtained. 
2 .   Parabolic  Interpolative  Eit  Approach 

A  more  recent  approach  (1975)  is  one  proposed  by 
Chang  et  al,  LRef.  6].   Chang  begins  by  assuming  a  surface 
current  on  rhe  antenna  (as  opposed  to  a  line  current 
approximation)  and  defines  a  new  Kernel  (Green's)  Function 
to  account  for  axial  variations  in  current.   He  also  derives 
a  different  expression  for  E^  from  the  relationship: 
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which  he  claims  is  more  accurate,  since  continuity  of  higher 
order  derivatives  of  'I'  do  not  have  to  be  specifically- 
considered  at  segment  boundaries  when  solving  the  Integral 
Equation.   He  then  assumes  a  Parabolic  Interpolative  Fit  for 
the  current  distribution  over  each  segment  (the  distribution 
over  the  whole  antenna  is  not  necessarily  parabolic)  where 
precise  values  of  current  are  exactly  defined  at  sample  points. 
He  also  defines  three  moment  functions  that  describe  the 
Kernel  Function  and  its  variation  with  respect  to  $  and  Z 
over  a  segment.   The  non-zero  field  points  are  then  expressed 
as  double  summations  of  the  moment  functions  over  each  inter- 
val.  To  avoid  the  double  summations,  he  splits  the  Kernel 
Function  into  two  factors  (a  technique  used  by  Hallen,  as 
well)  which  facilitates  the  computation  of  the  dominant  term 
(the  result  of  which  is  a  complete  elliptic  integral  of  the 
first  kind)  ,  then  redefines  the  Kernel  as  an  elliptic  integral 
plus  an  approximate  remainder.   With  this  Kernel  he  then 
computes  The  moment  functions  defined  above  while  redefining 
a  remainder  term  for  three  different  special  situations. 
These  moment  functions  can  now  Oe  evaluated  numerically 
through  adaptive  integration  techniques,  such  as  Gauss-Legendre 
Quadrature,  by  splitting  the  integrand  functions  into  Real  and 
Imaginarv  parts. 
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3 .   BesseL  Function  Fit  Approach 

Another  recent  approach  (19  81)  is  that  proposed  by 
Balzano  et  al ,  [Ref.  7],   Balzano  also  begins  by  redefining 
the  Kernel  Function  to  account  for  azimuthal  variations,  but 
casts  it  into  the  form  of  an  integral  of  a  Hankei  Function  of 
the  second  kind,  which  also  can  be  represented  as  an  infinite 
summation  of  Bessel  Functions.   The  magnetic  vector  potential 
is  then  expressed  as  a  double  integration  of  the  current 
density  on  the  conductor  times  the  integral  of  the  infinite 
summation  of  Bessel  Functions  described  above.   In  this  form, 
all  that  is  needed  is  an  accurate  description  of  the  axial 
current  distribution,  I(zr)  over  the  cylindrical  antenna. 

From  the  relationship  of  the  vector  potential  with 
the  E  and  H  fields,  an  expression  can  be  written  for  E_  as  a 
function  of  the  integral  of  I(z): 


E  (a,z)  =  -  Ar  f      J  <*'  -<2  -;2)H(2)  (a/  k2  -  ?2  ) 
z   '        8  irk       o         s    o  ^ 
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(k2  -  ;2)(  /   I(zr)e^Z  dz')e':<;Zd?      (23) 


I(z')  is  expressed  by  a  suitable  linear  combination  of 

orthonormal  functions  which  satisfy  the  boundary  conditions 

for  the  current.   Then  the  current  expressed  as  an  infinite 

sum  of  these  orthonormal  functions  is  equated  to  the  Fourier 

Transform  of  the  current  expressed  as  a  function  of  the 

vector  porential : 
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where:   e  (?  )  is  "the  Fourier  Transform  of  the  current 
n 

harmonics,  A   is  an  infinite  series  of  constants.   Eauaticn 
n 

(24)  is  substituted  into  Equation  (23)  and  solved  by  applying 

Galerkins  Method."   This  yields  an  infinite  system  of 

equations  which  can  be  truncated  and  solved  for  the  unknown 

A  .   A  similar  treatment  can  be  aoolied  to  the  other  non-zero 
n 

field  equations  to  cast  them  into  solvable  form.   However, 
the  equations  so  defined  converge  slowly;  therefore,  addi- 
tional manipulation  is  required.   The  additional  manipulation 
becomes  rather  involved  and  will  not  be  discussed  here. 
Details  can  be  found  in  Ref.  7. 
4 .   Discussion 

All  three  of  the  above  methods  are  theoretically 
sound,  all  three  give  reasonable  results  over  large  portions 
of  the  near-field  region,  and  ail  Three  lend  themselves  to 
numerical  solution  by  computer.   Of  concern  here  is  the 
amount  of  valuable  computer  resources  (time,  core,  and  cost) 
consumed  in  making  these  calculations.   Balzano  indicates  the 


*Galerkin' s  Method  solves  an  Integral  Equation  via  the  Method 
of  Moments  to  be  discussed  below.   Specifically,  Galerkin's 
Method  sets  the  Weighting  Functions  equal  to  the  3asis 
Functions  r 
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expressions  he  derives  for  field  components  are  slow  in 

converging,  and  even  after  modifications  to  his  basic 

derivation  to  hasten  the  convergence,  he  states  that  more 

than  1500  harmonics  (a  1500th  order  linear  system  of  equations) 

are  needed  to  satisfactorily  march  boundary  conditions. 

Chang  does  not  indicate  directly  the  amount  of 

computer  resources  consumed  by  his  method,  but  he  does  imply 

his  method  is  adopted  because  a  simpler  method  (The  Method  of 

Moments)  is  inaccurate: 

"...in  the  numerical  computation,  the  integral  from  -h  to  h 
is  usually  subdivided  into  a  finite  number  of  segments;  in 
each  segment  the  current  is  interpolated  by  a  polynomial 
with  coefficients  expressible  in  terms  of  values  of  current 
at  a  few  sample  points  within  that  segment.   The  higher 
derivatives  of  I(z)  are  therefore  discontinuous  at  boundary 
points  between  any  two  segments."  [Ref.  6] 

One  of  the  intentions  of  this  study  is  to  show  that 

the  Numerical  Electromagnetics  Code  (NEC),  which  uses  the 

Method  of  Moments,  has  been  modified  to  account  for  the 

discontinuities  of  the  first  derivative  at  segment  boundaries, 

and  that  this  modification  is  sufficient  to  yield  results  that 

are  accurate  for  engineering  analysis. 
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III.   DESCRIPTION  OF  NEC 

NEC,  Numerical  Electromagnetics  Code,  is  a  computer  pro- 
gram for  the  analysis  of  the  electromagnetic  response  of 
antennas  and  other  metal  structures.   It  computes  a  numerical 
solution  to  integral  equations  that  describe  the  currents 
induced  on  a  structure  by  voltage  or  current  generators  and/or 
incident  fields.   It  is  the  latest  in  a  series  of  programs. 
The  first,  called  3RACT ,  was  developed  by  MBAssociates  under 
funding  from  the  Air  Force  Space  and  Missile  Systems  Command. 
This  was  followed  by  AMP,  again  developed  by  MBAssociates, 
this  time  under  an  Office  of  Naval  Research  contract.   NEC 
evolved  from  AMP  and  was  developed  by  Lawrence  Livermore 
Laboratories  with  funding  from  the  Naval  Electronics  Systems 
Command.   It  is  a  user  friendly  cede  which  incorporates  the 
following: 

A.   FEATURES  OF  THE  CODE 

An  Integral  Equation  (I.E.)  for  modeling  the  current  on 
smooth  surfaces  Is  combined  with  an  I.E.  for  current  modeling 
on  thin  wires  to  describe  the  electromagnetic  response  of  an 
arbitrary  structure.   The  structure,  or  parts  of  it,  may  be 
active  or  passive,  located  over  a  perfect  or  imperfect  ground 
plane,  end  may  have  lumped-element  or  distributed  loading. 
Excitation  may  be  from  voltage  or  current  sources  on  the 
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structure,  an  incident  plane  wave  with  either  linear  or 
elliptic  polarization,  or  the  field  may  be  due  to  a  Hertzian 
dipole.   Output  may  include  current  and  charge  density,  power 
gain  or  directive  gain,  near-  or  far-zone  electric  or  magnetic 
fields,  impedance  or  admittance,  total  radiated  power  or  input 
power.   It  is  suitable  for  either  antenna  analysis  or  cross- 
section  scattering  and  electromagnetic  pulse  (EMP)  studies. 
The  code  will  handle  junctions  of  wires  of  uneven  radius,  free 
ends  of  wires  that  have  finite  radius,  wires  of  variable  radius, 
and  coupling  between  wires. 

NEC  utilizes  the  Gauss-Doolittle  Method  for  solving  the 
matrix  equations  generated  by  the  Method  of  Moments  when 
solving  the  integral  equations.   It  also  allows  for  use  of 
rotational  or  plane  symmetry  to  reduce  computation  time.   When 
the  impedance  matrix  is  too  large  to  be  contained  in  core ,  NEC 
has  the  option  of  storing  portions  out  of  core.   It  allows  the 
1  self -interaction  matrix'  for  a  structure  to  be  computed, 
factored  for  solution,  and  stored  on  tape  or  file.   A  solution 
for  a  new  antenna  that  enters  this  environment  requires  only 
the  evaluation  of  the  '  self -interaction  matrix'  for  the 
antenna,  the  mutual  antenna  to  environment  interactions,  and 
matrix  manipulations  for  a  partitional  matrix  solution. 

B.   ZONING  CONSIDERATIONS 

NEC  is  a  discrete  sampling  code  where  a  complex  structure 
must  be  dissected  into  a  number  of  simole  elements  (wires  or 
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plates)  to  which  the  Electric  Field  Integral  Equation  (EFIE) 
or  Magnetic  Field  Integral  Equation  (MFIE)  are  applied.   As 
with  any  description  of  the  real  world,  there  are  approxima- 
tions, but  because  of  the  versatility  in  modeling  the  geo- 
metry of  a  structure,  its  approximations  more  closely 
resemble  nature.   The  resemblance  is  strongly  influenced  by 
the  choice  of  zoning  (i.e.  dissecting)  the  structure  in  the 
program.   The  smaller  the  geometric  elements,  the  closer  the 
model  comes  to  reality.   However,  the  smaller  the  elements, 
the  larger  the  number  of  elements ,  which  means  the  larger  the 
matrix  of  equations  and  hence,  the  more  costly  the  solution. 
There  is  a  point  beyond  which  smaller  zones  do  not  yield  a 
substantially  more  accurate  solution;  it  may  be  necessary  to 
rerun  the  problem  with  increasingly  smaller  elements  to  find 
the  point  of  diminishing  returns.   The  choice  of  proper  zoning 
then  is  gained  with  experience  and  becomes  as  much  of  an  art 
as  it  is  a  science.   The  guidelines  for  the  science  asoect 

as  follow: 

I .   Wires 

Segments  should  follow  the  paths  of  conductors  using 
a  piece-wise  linear  fit  on  curves.   Generally,  segment  lengths 
(a)  should  be  less  than  .IX;  shorter  segments  (.0  5\)  or  less 

may  be  needed  at  critical  regions  (junctions  or  curves). 

_3 
Segments  smaller  than  10  \    should  be  avoided  since  the 

similarity  of  constant  and  cosine  components  lead  to 

numerical  inaccuracy.   The  radius  of  the  wire  (a)  relative  to 


A  depends  on  the  Kernel  used  in  the  Integral  Equation.   Two 
options  exist.   The  thin-wire  Kernel  models  a  filament  current, 
while  the  extended  Kernel  models  a  uniform  current  distribu- 
tion around  the  segment  surface.   The  field  of  the  distribu- 
ted current  is  approximated  by  the  first  two  terms  in  a 
series  expansion  of  the  exact  field,  in  powers  of  a".   The 

n 

first  (a  )  term  is  identical  to  the  thin-wire  Kernel;  the 
second  term  extends  the  accuracy  for  larger  values  of  a. 
Both  Kernels  incorporate  the  thin-wire  approximations  (see 
Section  II. A)  and  both  require  2ira/X<<l.   The  thin-wire 
Kernel  requires  a  A/a>3;  the  expended  Kernel  relaxes  this 
to  A/a>2.   These  values  ensure  errors  are  less  than  1%.   The 
extended  Kernel  is  used  at  free  wire  ends  and  between  parallel, 
connected  segments.   The  thin-wire  Kernel  is  always  used  at 
bends . 

2 .   Surfaces 

A  conducting  surface  is  modeled  by  small  flat  surface 
patches  which  conform  as  closely  as  possible  to  curved  surfaces 
The  parameter  defining  a  patch  is  a  normal  unit  vector, 
originating  from  the  center  of  the  patch,  defined  in  Cartesian 
coordinates.   Each  patch  must  be  connected  by  a  wire  at  the 
patch  center  for  the  program  to  integrate  the  surface  current. 
The  code  divides  each  patch  into  four  equal  patches  about  the 
wire  end,  along  the  unit  vector  lines  describing  The  surface 
of  the  patch.   An  interpolation  function  is    applied  to  the 
four  patches  to  represent  placement  of  the  current  onto  the 
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patches,  and  the  function  is  numerically  integrated.   Patches 
with  wires  connected  to  them  (active  patches)  should  be 
chosen  to  be  approximately  square  with  sides  parallel  to  the 
unit  vectors  defining  the  surface.   Only  one  wire  may  connect 
to  a  patch,  that  wire  may  not  be  connected  to  another  patch, 
and  it  may  not  lie  in  the  plane  of  the  patch.   A  minimum  of 
about  2  5  patches  should  be  used,  per  wavelength  of  surface 
area;  the  maximum  size  for  an  individual  patch  is  about  .  GLi- 
square  wavelengths.   The  number  of  patches  used  increases, 
and  the  size  of  each  patch  decreases,  as  the  radius  of 
curvature  decreases.   Smaller  patches  should  be  used  at 
edges  since  the  current  amplitude  may  vary  rapidly  in  this 
region.   Long  narrow  patches  should  be  avoided.   Patches  are 
restricted  to  modeling  voluminous  bodies  with  closed  surfaces; 
parallel  surfaces  on  opposite  sides  cannot  be  too  close 
together . 

3 .   Ground  Plane 

For  a  perfectly  conducting  ground,  the  code    generates 
a  reflected  image.   Structures  may  be  close  to,  or  contact, 


the  ground;  however,  for  a  horizontal  wire:   /  h   +  a  >10   \ 
where  a  =  wire  radius ,  h  -  height  of  wire  axis  above  the 
ground  plane,  and  h/a>3. 

A  finitely  conducting  ground  may  be  modeled  by  an 
image  modified  by  the  Fresnei  Plane-wave  Reflection 
coefficients.   This  method  is  fast,  but  of  limited  accuracy 
and  should  not  be  used  for  structures  close  to  the  ground,  or 
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having  a  large  horizontal  extent  over  the  ground.   The 
Sommerf eld/Norton  model  uses  the  exact  solution  and  is 
accurate  close  to  the  ground;  the  horizontal  restriction  is 
the  same  as  for  a  perfect  ground.   This  model  is  only  used 
for  wire-wire  interactions,  for  surfaces  the  code  reverts  to 
Fresnel  Reflection  coefficients.   At  the  present  time,  a 
ground  stake  cannot  be  used,  but  wires  may  end  on  a  perfectly 
conducting  ground  plane  if  the  derivative  of  the  current  at 
the  end  of  zhe   wire  is  zero.   For  wires,  there  are  options  for 
a  radial-wire  ground  screen  approximation  and  a  two-medium 
ground  approximation  based  on  modified  reflection  coefficients 
These  options  allow  considerable  savings  in  computational  time 
when  limited  accuracy  can  be  tolerated. 

C.   INTEGRAL  EQUATIONS 

1.   Electric  Field  Integral  Ecuation  (EFIE) 


The  EFIE  is  used  for  thin-wire  structures  with  small 
conductor  volume.   It  is  derived  from  the  electric  field 
representation  for  a  current  distribution  confined  to  the 
surface  of  a  perfectly  conducting  body: 


E(p)  =  _n  JL-   $     J    (r?)  G(r.r')  dA'  (25) 

-'     ^TC   s   s 


where:   t'  is  the  source  point,  r  is  the  observation  point, 

2= 


s  a  surface  current  density,  G(r,r')  =  ( 3  I  +■  VV)g(r,: 


t  > 


_l Q I p-p '     I        i  .       .. 

(r,rT)    =    e    JhJ]  '/!r-rTj  ,    3    =    phase   velocity,    n    =    intrinsic 
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impedance  of  medium.   The  principal  value  integral  /   is 
indicated  so  r-*r?  in  the  limit. 

An  integral  equation  for  the  current  induced  on  S  by 
an  incident  field  E   can  be  obtained  from  Equation  (25)  and 
the  boundary  condition: 


n(r)  X[E~  (r)  +  Es  (r)]  =  0  (26) 


where  rz S ,  n(r)  is  the  unit  normal  vector  at  r,  Eb  is  the 

field  due  to  the  induced  current  J  .   This  integral  equation 

s  &      m 

is  given  by: 


-n(r)  X  EX(r)  =  -j  JL-  n(r)  X  J   J  (rf) 

M-ttB         *s  s 


(62I  +  77)g(r,r? )  dA'  (2?) 


This  vector  integral  equation  is  reduced  to  a  scalar  integral 
equation  for  the  thin-wire  model.   With  the  surface  current 
density  replaced  by  a  filament  current  and  the  boundary 
condition  enforced  in  the  axial  direction,  this  gives: 


.2 

-S-E^CF)  =  -Jtt~  J"  KsT)(B2s.sT  -  ,3,  ,)g(F,F,)ds'   (28) 
^•n-gj  3  s9  s  ' 
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where:   s  =  distance  along  the  wire  axis  at  r,  and  s'  =  unit 
vector  along  the  wire  axis.   This  is  the  EFIE  solved  by  NEC. 
2 .   Magnetic  Field  Integral  Equation  (MF1E) 


The  MFIE  is  derived  from  the  integral  representation 
for  the  Magnetic  field  of  a  surface  current  distribution  J  : 


Hs(r)  =  ~-  f   J   (r)  X  t-^-t-   g(r,r')  dAf  (29) 

M-TT      s        dr 
5 


If  J   is  induced  bv  an  external  incident  field  H  ,  then  the 
s 

total  magnetic  field  inside  the  perfectly  conducting  surface 
must  be: 


H1(r)  +  HS(r)  =  0  (3G) 


where:   H   is  the  incident  field,  and  H   is  the  scattered 

field.   Let  r->r   along  n(r  )  from  inside  the  surface,  then 

o      °  o 

Equation  (30)  with  Equation  (29)  substituted  becomes: 


•n'r  )  X  Hx(r  )  =  nCr  )  X  i—   lim   /  J  (r') 

O  O  O       HIT  S 

r-»-r  s 
o 


X  -^-r  g  (  r ,  r  *  )  dA '  (31) 

3r '  ^ 


The  limit  is  evaluated  using  a  result  of  potential  theory  to 
yield: 


uo 


■n(r  )  X  H1(r  ) 
o        o 


-1/2 J  (r  )  +  f-  i      n(F  ) 

SO       M-TT    S       O 


X[J  (r1)  X  -2-,-  g(r  ,r'  )]  dAf 
s        3r'  6   o' 


(32) 


This  vector  is  resolved  into  two  scalar  equations  along  the 
orthogonal  surface  vectors  t   and  t  ,  where: 


t, (r  )  X  t„(r  )  =  n(r  ) 

1   o      2   o        G 


(33) 


which  are  given  by 


t0(r  )-H1(r  )  = 
2   o      o 


-1/2  t,  (r  )  -J  (r  ) 
1   o    so 


J7-   /      t0(r    )-CJ    (rf)    X   J-r  g(r    ,r')]dA' 
Utts2o  s  3rT3o 


(34) 


and 


■t,  (r    )  •H1(r    )    = 
1      o  o 


•1/2    t0(r    )  -J    (r    )    + 
2      o         so 


+   i-  $      t,  (r    )'[J    (r')    X   -$-r  g(r    ,r')]dAr 
M-ttsIo  s  3  r  o 


(35) 
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3.   Hybrid  EFIE/MTIE 

NEC  models  thin  wires  using  EFIE  and  for  surfaces 

uses  MFIE.   For  a  structure  with  both  wires  and  surfaces,  r 

in  Equation  (28)  is  restricted  to  the  wires,  with  the  integral 

for  Es(r)  extending  over  the  complete  structure.   Hence, 

Equation  (28)  is  used  for  wires,  while  the  more  general  form, 

Equation  (27)  is  used  for  surfaces.   Likewise,  r   is 

o 

restricted  to  surfaces  in  Equation's  (3  4)  and  (35),  with  the 

p  

integrals  for  H  (r)  extending  over  the  complete  structure. 

The  resulting  coupled  integral  equations  are: 
For  r  on  wire  surfaces : 

-s-E^r)    =    -j    — ^     /    I(s?  ) -Cs2s-s?    -      ^       t]g(r,r')    ds  ? 

471dt  3S3S 

Li 


•j    T—r        /      J    (r)-[32s-r-^-r    •    |~lg(r,r')    dA» 
J    M-ttS  s  3r '  8s    a 
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(36) 


For   r   on    surfaces    excluding   wires 


t„(r)    HX(r)    =    -  i-  t9(r)-    /    Ks'JCs'    X  -ij-  g(r,rT)]    ds 
2  M-tt       2  T  3rT 

Li 


1/2    t,  (rW    (r)    -   ~   /        t«(r) 


i_  ,,       2    --> 

■IT 


"1 


[J    (rf)    X  -2-r  gCr,r')]    dAf  (37) 
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and 


-t,  (r)  -Hx(r)    =    £-  tn  (r)     •    /    I(s')[s?    X   ^  g(r,r?)]    ds ' 

1  4tt       1  .  3r 


■l/2t0(rW    (r)    +   i-   r        t.  (r) 

2  3  41         S,  1 


[J    (rf)    X  r-r  g(r,r»)]   dAf  (38) 


D.   METHOD  OF  MOMENTS  AND  THE  MATCH  POINT 

The  Method  of  Moments  is  a  technique  whereby  an  integral 
equation  is  reduced  to  a  system  of  linear  algebraic  equations 
which  are  easily  handled  by  high  speed  digital  computers. 

1 .   Mathematical  Concept 

The  method  applies  to  an  inhomogeneous  linear 
operator  of  the  form: 


(39) 


where  L  is  a  linear  operator  with  domain  DT  containing  the 
function  f,  which  is  an  unknown  response  to  an  excitation  e, 
which  is  given,  and  found  in  the  range  of  1. 

If  a  solution  to  Equation  (39)  exists  and  is  unique 
for  all  e,  then  the  inverse  operator,  L   ,  exists  such  that: 

r  =  L"1?  (uo) 

which  represents  a  solution  to  Equation  (39). 
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The  unknown  function,  f,  may  be  expanded  in  a  series  of 
N  linearly  independent  Basis  Functions  spanning  DT  as: 

N 

f  =  Z      a.  F.  (^1) 

Substituting  Equation  (41)  into  Equation  (39),  and  taking 
the  inner  product  with  a  set  of  linearly  independent 
weighting  functions  Cw.},  defined  in  the  range  L  and  spanning 
sub-space  Sr  ,  results  in  a  set  of  equations  for  the 
coefficients,  a.,  of  Equation  (41).   This  set  of  equations: 

N 

Z      a.<w.,  Lf.>  =  <w.,  e>     i  =  1,  2,...N  (42) 

.-,111       i 
j  =1   J        J 

can  be  written  in  matrix  notation  as: 

EG]  CA]  =  [E]  (43) 


where:   G..  =  <w.  ,  Lf . >  ,  A.  =  a.,  and  E.  =  <w. -  e> .   This 
i]      i'    3     3     3        l      i' 

procedure  projects  the  space  of  e  onto  the  sub-space  Sr  ,  and 
enforces  the  equality  between  projections. 

Hence,  if  a  solution  to  Equation  (43)  exists  and  is 

_  _  i 

unique  for  all  e,  then  the  inverse  operator,  L   ,  also 

exists  such  than: 
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[A]  =  [G]"1  [E]  (UU) 


which  is  a  solution  to  Equation  (43). 

The  efficiency  of  computations  and  accuracy  cf 
solution  is  largely  dependent  on  the  choice  of  Basis  Function. 
Factors  which  should  guide  this  choice  are: 

a)  accuracy  of  solution  desired 

b)  ease  of  evaluation  of  matrix  elements 

c)  size  of  the  matrix  that  can  be  inverted 

d)  the  realization  of  a  'well-conditioned'  matrix. 
There  are  two  general  classes  of  Basis  Functions  from  which 

to  choose,  entire  domain  and  sub-domain.   The  sub-domain  class 
has  fewer  elements;  hence,  its  execution  time  is  less,  it 
simplifies  the  evaluation  of  the  inner  product  integral  and 
still  ensures  the  matrix  [Q]  will  be  'well-conditioned'. 
2  .   NEC's  Application 

NEC's  interpolation  of  current  on  wires  enforces 
continuity  of  current  and  linear  charge  density  at  segment 
ends.   This  imposes  additional  restrictions  on  the  current 
and  its  derivative  with  respect  to  distance  along  the 
conductor.   NEC  forces  the  current  on  the  conductor  to  match 
the  integral  equation  at  specific  points;  these  'match  points' 
are  located  at  the  center  of  each  segment. 

NEC's  application  cf  the  Method  of  Moments  starts 
with  the  Integral  Equations:   Equation's  (23)  and  (34)  to 
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(35),  or  Equation's  (36)  to  (38);  the  equation  selected 
depending  on  the  structure  being  modeled.   These  integral 
equations  are  the  basis  from  which  the  system  of  linear 
equations  will  be  generated. 

A  source  model  with  a  known  excitation  voltage  can 
be  related  to  the  E-field  impressed  on  the  antenna;  hence, 
E   (=  e  from  Equation  (39))  is  known.   We  know  we  are  locking 
for  a  more  accurate  description  of  the  current  distribution 
on  the  antenna;  hence,  the  current  I  (=  f  from  Equation  (39)) 
is  the  unknown  response.   The  operator  L  (an  integral 
operator  in  this  case)  is  the  composite  of  all  other  terms/ 
factors  (other  than  E  and  I)  found  in  the  integral  equations. 
Hence,  these  terms  can  be  cast  into  the  form  of  Equation  (39) 
as  : 

E1  =  L  I  (45) 

To  generate  the  system  linear  equations,  we  have 
divided  the  structure  to  be  modeled  into  a  number  of  segments 
according  to  the  zoning  guidelines  discussed  above.   The 
current  distribution  on  the  antenna  is  now  approximated  as 
the  summation  of  currents  found  en  each  segment.   That  is: 


M 
I  =  Z         I.  g.  C^+6  ) 

1=1    ]   ^ 
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which  has  the  same  form  as  Equation  (4-1).  Here  e.  are  the 
3asis  Functions  and  I.  is  the  system  of  coefficients.  For 
wires : 


I.g.  =  A.  +  B.  sin(B(s-s.)  +  C.  cos( 3 ( s-s . ) )*        (47) 
3  1     ]     ]  ]      ]  3 


where:   s.  =  center  of  the  ith  sub-section,  and  3  =  free 
space  wave  number.   As  a  result  of  current  continuity 
considerations  at  segment  ends,  each  sub-section  will  have 
only  one  unknown. 

For  patches,  the  surface  current  density  is: 


J  (r)  =  u  K  .  +  v  K  .  (43) 

s         u]       v- 


where:   K  .  and  K  .  are  vector  components  of  the  surface 
current  density,  J  (r),  on  each  patch.   Each  K.  is  constant 
over  the  patch,  which  is  usually  rectangular,  but  can  be  any 
reasonable  shape.   Again,  current  continuity  considerations 
for  the  patch  yield  only  two  unknowns,  one  for  each  direction 

3ecause  of  the  nature  of  the  Kernel  functions  in  the 
integral  equations,  the  choice  of  basis  functions  is  mere 
critical  on  wires  than  on  surfaces. 


"This  form  has  been  found  to  provide  rapid  solution  conver- 
gence and  fields  of  sinusoidal  currents  are  easily  evaluated 
in  closed  form. 


For  simplicity,  the  rest  of  this  section  will  be 
restricted  to  the  case  for  only  wires,  as  this  is  in 
consonance  with  the  system  for  this  study.   Assuming  the 
wire  is  aligned  with  the  z-axis  (as  it  is  in  this  study), 
Equation  (45)  can  be  re-written  as: 

E1(z)  =  L  I  (z')  (49) 

where  the  zf  indicates  the  current  is  on  the  surface  of  the 
wire.   Equation  (46)  can  be  rewritten  as: 


N 
I(zf)=Zl.  g.(z')  (50) 

5=1  ^  ^ 


Equation  (50)  defines  the  set  of  sub-domain  basis  functions 

to  be  used  in  computations.   For  finite  M,  the  sum  of  I.  g.(zf) 

cannot  exactly  equal  the  general  current  distribution*",  so 

the  I.  g.(z')  need  to  be  chosen  as  close  to  the  actual 

]   D 

distribution  as  possible.   This  is  accomplished  by  restricting 
the  suDtort  for  I.  g.(z')  to  a  localized  region  of  the  surface 
near  the  center  of  the  segment.   Each  basis  function  has  a 
peak  on  the  segment  on  which  it  is  defined  and  extends  only 


"The  approximation  can  be  brought  arbitrarily   close  to  the 
exact  solution  by  choosing  N  sufficiently  large.   Fortunately, 
reasonable  accuracy  can  be  obtained,  in  most  cases,  with  a 
relatively  small  number  of  segments  (N  -  10/X). 
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onto  segments  to  which  It  is  directly  connected,  going  to 
zero  with  zero  derivative  at  the  outer  ends  of  connected 
segments . 

Equation  (50)  is  then  substituted  into  Equation  (M-9) 
to  give: 


M  N 

E1(z)  =  LL  I   I   g.(z')]  =   S   I.  LCg.(z')]  (51) 

j=l   J   3        j=1  1  1 


where:   I.  and  the  summation  have  been  moved  outside  the 
] 

domain  of  L  as  a  result  of  its  linearity. 

The  next  task  is  to  identify  a  set  of  weighting 

functions,  (w.).   When  the  weighting  functions  are  set  equal 

to  the  basis  functions  (i.e.  w.  =  g.),  this  is  known  as 

i     l 

Galerkins  Method."   In  NEC,  w.  is  chosen  as  a  set  of  Dirac 

i 

Delta  Functions: 


w.(z)  =  6(z-  *i )  (52) 

l  2 


A  1 
where:   y—  is  the  set  of  match  points  found  at  the  center  o: 

each  segment  Ai.   Use  of  this  weighting  function  results  in 

a  point  sampling  of  the  integral  equations  known  as  the 

"collocation  method"  of  solution. 


'Employed  in  Hall en's  approach  (See  Section  II.C.l) 


^9 


Using  this  weighting  function  (Equation  (52))  and 
taking  the  inner  product  of  each  side  of  Equation  (51)  gives 


i        N 
<w.,E(z)>=  z      I.<w.,L[g(z')]>  (5  3) 


which  has  the  form: 

[V]  =  [I]  [Z]  (54) 

where:   [V]  is  a  matrix  of  weighting  functions  applied  to  the 
known  impressed  field,  [I]  is  a  matrix  of  unknown  coefficients, 
and  [Z]  is  a  mapping  matrix  of  weighting  functions  applied  to 
the  assumed  basis  functions. 

From  previous  arguments,  we  know  that  if  a  solution 
to  Equation  (54)  exists  and  is  unique  for  all  E  ,  then 
Equation  (54-)  can  be  transposed  to: 


[I]  =  [Z]"1  [V]  (55) 


In  this  form,  the  system  of  equations  can  be  solved  by 
imposing  the  boundary  condition: 


i1  +  Es  =  0  '56) 

z     z 
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on  E  ,  which  then  states  that  the  tangential  scattered  field, 

E   (which  generates  the  current) ,  is  equal  to  the  negative 

of  the  known  impressed  field,  E  .   Stated  in  eauaticn  form: 

'   z 


E1  =  -ES  (57) 

z      z 


Hence,  the  system  of  equations  can  be  solved  for  the  exact 
current  at  the  march  points,  and  by  enforcing  current  and 
charge  continuity  at  segment  ends  ,  a  very  close  approximation 
to  the  actual  current  distribution  can  be  obtained. 

Equation  (54)  is  solved  by  the  Gauss-Doolittle  method. 
In  this  method  the  matrix  [Z]  is  defined  as  two  sub-matrices: 


[Z]  =  Lp  (59) 


where  sub-matrix  [A]  with  elements  A.,  represents  the  electric 
field  at  the  center  of  segment  i  due  to  the  jth  basis 
function,  centered  on  segment  j.   Sub-matrix  L3] ,  with  elements 
B..,  represents  magnetic  fields  generated  by  segment  basis 

-'-J 

functions.   In  the  execution  of  this  method,  the  matrix  [Z] 
is  factored  into  the  product  of  an  upper  triangular  matrix 
[U]  and  a  lower  triangular  matrix  [L],   That  is: 

[Z]  =  CL]  [U]  (59) 
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E a nation  (54)  then  becomes 


[V]  =  [I]  CL]  CU]  (60) 


and  [I]  is  comDuted  in  two  stages 


CV]  =  [L]  [F]  (51) 


and 


CF]  =  CU]  EH  (52) 


where  Equation  (61)  is  solved  first  for  [F]  (an  intermediate 

matrix)  by  forward  substitution,  and  then  for  [I]  by  backward 

substitution. 

Now  that  the  current  has  been  determined,  the  other 

non-zero  scalar  field  components  (E   and  Hx)  can  easily  be 

r       $ 

computed . 
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IV.   MODEL  DESCRIPTION  AND  COMPUTATIONAL  RESULTS 

A.   MODEL  DESCRIPTION 

The  basic  model  used  for  this  study  is  a  simDle  monooole 
colinear  with  the  Z-axis  and  mounted  perpendicular  to  a  per- 
fectly reflecting  ground  plane,  that  is  itself  located  in  the 
X-Y  plane.   As  can  be  seen  from  Figure  1,  the  antenna  is  7  5 
meters  in  height  and  .5  meters  in  diameter.   It  is  radiating 
at  a  frequency  of  1  Megahertz  (  the  middle  of  the  AM  broad- 
cast band)  with  a  wavelength  (X=c/f)  of  3  00  meters;  it  is  a 
quarter-wavelength  monopole  operating  at  its  resonant  fre- 
quency.  Transmit  power  is  fed  to  the  model  through  its  base 
by  a  voltage  that  is  normalized  to  obtain  10 CO  watts  of 
radiated  power. 

Virtually  all  modeling  schemes  start  from  the  basic 
premise  that  E-  and  H-field  component  values  can  be  determined 
from  a  knowledge  of  the  current  distribution  on  the  antenna. 
Three  generalized  methods  for  describing  the  current  distri- 
bution will  be  discussed  here:   1)  assuming  a  current 
distribution,  2)  modal  modeling,  and  3)  sub-domain  mcdeling. 

1 .   Assumed  Current  Distribution 

This  is  the  classical  method  as  described  in  [Ref.  M-j. 
Typically,  the  assumed  distribution  is  sinusoidal,  which 
results  in  a  tangential  E-field  component  substantially 
different  from  zero.   This  may  not  cause  a  significant  error  in 
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the  far-field ,  but  in  the  near-field  where  stored  energy 
strongly  depends  on  the  current  shaoe,  large  errors 
generally  occur. 

2 .  Modal  Modeling 

To  improve  on  the  assumed  distribution,  the  modal 
modeling  method  was  devised.   By  this  method  the  current 
distribution  can  be  described  as  the  sum  of  an  infinite 
system  of  equations  which  precisely  describe  the  current 
distribution  in  the  limit.   This  is  the  approach  used  by 
Balzano  [Ref.  7].   Using  an  infinite  series  of  Bessel-Fourier 
Integrals,  which  can  be  truncated  to  allow  for  desired 
accuracy,  a  system  of  closed  form  equations  can  be  generated 
which  lend  themselves  to  computer  solution.   This  method 
enforces  boundary  conditions  (i.e.  E_  =  0)  over  the  surface 
of  the  antenna  with  precision  (in  the  limit),  but  there  is  a 
practical  limitation  that  must  be  imposed  on  the  size  of  the 
system  of  equations  that  precludes  an  exact  solution.   If  the 
current  is  expressed  by  a  sum  of  harmonics,  the  total  number 
of  harmonics  which  must  be  used  in  the  computations  to  match 
boundary  conditions  Is  related  to  the  ratio  of  dipole  height 
to  gap  size. 

3 .  Sub-domain   Modeling 

To  avoid  much  of  the  manhematicai  complexity  of  modal 
modeling,  sub-domain  modeling  divides  a  given  structure  into 
a  concaienaied  series  of  sub-structures,  or  elements.   The 
actual  current  distribuTicn  on  the  given  structure  is 
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approximated  by  geometric  distributions  (sin,  cos,  and  linear 
constants)  on  the  elements.   This  is  the  method  used  by  NEC 
[Ref.  2],   The  NEC  code  uses  the  scalar  form  of  the  integral 
equations  to  determine  the  current  in  each  of  the  elements, 
then  combines  the  results  from  each  substructure  to  deduce  the 
performance  of  the  total  structure. 

B.   MODEL  CONFIGURATIONS 

In  this  study  three  configurations  of  the  model  were 
assumed  and  values  for  the  tangential,  radial,  and  peak  (time 
and  spatial)  fields  were  calculated  for  each.   The  three 
configurations  were  for  a  structure  modeled  with  5,  15,  and 
2  5  sub-domains.   In  each  of  the  configurations,  a  null  value 
for  the  tangential  E-field  was  imposed  at  the  center  point  of 
each  sub-domain  (this  point  is  called  the  match  point).   Since 
a  closed  form  solvable  function  has  not  been  identified  which 
describes  the  current  on  the  antenna,  while  maintaining  a 
tangential  E-field  equal  to  zero,  it  has  become  necessary  to 
approximate  the  current  distribution.   This  is  done  in  NEC 
through  the  Method  of  Moments  as  discussed  in  Section  III.D. 
Using  this  method  the  current  can  be  approximated  to  an 
arbitrary  degree  of  accuracy  by  increasing  N.   Fortunately, 
reasonable  accuracy  is  obtained  in  most  situations  with  a 
relatively  small  number  of  segments,  typically  M  =  10/A. 
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C.   EXPECTATIONS/OBSERVATIONS/ANOMALIES 

This  study  utilized  three  configurations  to  determine  two 
relevant  facts:   a)  whether  the  calculated  value  of  the 
near-field  was  dependent  on  the  number  of  segments  used, 
and  b)  whether  the  size  of  the  segment  conformed  to  the 
guidelines  used  in  the  far-field. 

A  non-analytic  overview  of  the  pIots  generated  by  this 
study  indicated  that  the  calculated  value  of  the  near-field 
was  independent  of  the  number  of  segments  into  which  the 
structure  was  divided.   This  was  determined  by  overlaying 
plots  of  five  and  fifteen  segments  on  a  twenty-five  segment 
plot.   An  example  of  one  of  these  plots  is  shown  in  Figure  4. 
The  spurious  field  emitted  from  the  base  of  the  anrenna  has 
significantly  different  values.   These  variations  were  the 
result  of  using  different  values  of  input  voltage,  when  The 
configuration  was  changed,  to  maintain  a  normalized  radiated 
power  output  of  10  0  0  watts. *   The  values  along  The  antenna 
and  at  its  top  are  essentially  the  same,  for  respective 
positions  on  the  antenna,  for  all  three  configurations. 

The  NEC  manual  [Ref.  2]  recommends  segment  lengThs  (a) 
less  Than  X/10  for  normal  modeling,  and  for  critical  regions 
a  A  <  X/20  is  recommended.   For  a  wavelength  of  3  00  meters 
as  used  in  this  study,  segment  lengths  of  i  <  300/10  =  30 


"Radiated  power  depended  on  the  magnitude  of  the  E-Field  in 
each  segment,  which  in  Turn  depended  on  The  segment  length. 
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meters,  or  A  <  300/20  =  15  meters  are  suggested.   Lengths  used 
in  each  configuration  were:   A,  =  75/5  =  15  meters,  A~  =  7  5/15 
=  5  meters,  and  a  ~  =  75/25  =  3  meters.   Again,  from  a  non- 
analytic  overview  of  plots  similar  to  Figure  4,  it  was 
determined  that  reasonable  additional  resolution  was  obtained 
in  going  from  a  five  segment  to  a  fifteen  segment  configura- 
tion (i.e.  A  =  15  meters  to  A  =  5  meters),  but  the  increased 
resolution  in  going  from  the  fifteen  segment  to  the  twenty- 
five  segment  configuration  (i.e.  A  =  5  meters  to  A  =  3  meters) 
was  not  justified  by  the  increased  computer  resources  it 
consumed.   From  this  observation  it  could  be  concluded  that 
smaller  segment  lengths  are  required  for  near-field  calcu- 
lations, and  that  the  lengths  are  on  the  order  of  1/5  the 
value  specified  for  the  far-fields  (i.e.  A<X/50  for  normal 
situations  and  A<X/100  for  critical  areas).   Additional 
tests  and  models  should  be  run  to  modify,  or  confirm  this 
observation. 

Intensity  of  the  near-field  is  independent  of  the 
segment  size,  but  smaller  segment  size  is  required  to  get  a 
resolution  of  how  that  intensity  varied  as  a  function  of 
position  along  the  antenna.   Based  on  these  observations,  a 
determination  was  made  to  restrict  all  further  discussion 
to  the  fifteen  segment  configuration  to  avoid  duplication 
of  effort. 

In  the  course  of  the  investigation,  plots  of  the  E-  and 
H-fields  as  a  function  of  distance  along  the  Z-axis  were  run 
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from  the  base  of  the  monopoie  to  its  height,  and  also  from 
the  base  to  twice  the  height  of  the  monopoie.   The  plots  to 
just  the  height  of  the  monopoie  were  run  to  obtain  increased 
resolution  as  a  function  of  position  along  the  monopoie. 
Examination  of  the  single  and  double  height  plots,  side  by 
side,  revealed  no  increased  resolution  or  additional  infor- 
mation was  obtained  from  a  single  height  olot  that  was  not 
available  on  the  double  height  plot,  so  only  plots  to  twice 
the  antenna  height  are  discussed  in  this  section. 

Additionally,  closely  spaced  plots  within  a  segment  were 
run  to  determine  if  there  was  any  variation  as  a  function  of 
intra-segment  position.   This  type  of  plot  was  run  at  the 
center  of  the  antenna  for:   a)  the  end  of  a  segment,  b)  the 
center  of  a  segment,  and  c)  a  position  half-way  between  the 
positions  described  in  a)  and  b).   This  type  of  plot  was  also 
run  for:   a)  the  top  of  the  antenna,  b)  a  quarter  of  a  segment 
above  and  below  the  top,  and  c)  a  half  a  segment  above  and 
below  the  too.   In  both  cases,  there  was  no  observable 
difference  in  the  a.'),    b)  ,  or  c)  plots.   Consequently,  they 
have  not  been  included  in  the  discussion  of  plots  that 
follows . 

As  a  result  of  the  above  culling,  the  following  are  the 

plots  included  in  Appendix  A: 

I.   Z-field  olots  (E  ,  E  ,  E    .  )  as  a  function  of  Z  to 

2    r    peak 

twice  the  height  of  the  monopoie,  at  distances  of    1,  2,  5, 
10,  100,  and  1000  radii  from  the  axis. 
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2.  E-Field  plots  (E  ,  E  ,  E    ,  )  as  a  function  or  R  to 

z    r   peak 

a  distance  of  one  (1)  wavelength  from  The  monopole ,  for 
heights  of  2/100ths,  1/4,  1/2,  1,  and  2  rimes  the  height 
of  the  monopole. 

3.  H-field  plots  (H)  as  a  function  of  Z  to  twice  the 
height  of  the  monopole,  at  distances  of  1,  2,  5,  13,  100, 
and  10  0  0  radii  from  the  axis. 

4.  H-field  plots  (H)  as  a  function  of  R  to  a  distance  of 
one  (1)  wavelength  for  heights  of  2/100ths,  1/4,  1/2,  1, 
and  2  times  The  height  of  the  monopole. 

Each  of  these  categories  are  discussed  in  further  detail 
below. 

1.   E(Z)  Plots 

As  shown  in  Figure  5,  Z^(Z)  along  the  antenna  surface 
was  much  lower  in  magnitude  (about  6  0  dB  lower)  than  its  value 
at  the  too.   These  lower  values  varied  bv  an  order  of  magni- 
tude  as  a  result  of  NEC  forcing  E   =  0  at  match  points  and 
approximating  values  (from  the  current  distribution)  else- 
where.  These  radical  variations  were  quickly  attenuated  ar 
the  observation  point  (O.P.)  moved  awa2/  from  the  surface.   At 
a  distance  of  five  radii  (Figure  5),  the  integrating  effect 
of  a  wider  field  of  view  had  essentially  smoothed  cut  the 
curve,  and  also  decreased  the  difference  in  magnitude  between 
values  at  the  top  of  the  antenna  and  Those  along  its  length. 
When  the  O.P.  moves  to  a  distance  of  10  0  radii,  the  integrating 
effect  will  have  created  an  almost  uniform  magnitude  along 
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Plot  Ez  vi  Z  for  a  model  of  15  togmontt. 
Plot  It  ot  on*  radlut.  tho  Surfoco  of  tho  Monopoly. 
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Plot  Ez  v«  Z  for  a  modol  of  15  togmonfs. 
Plot  It  ot  t«n  radii.  2.7  M«tor«  obevo  »ho  Surfoco. 
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the  antenna,  and  this  uniform  magnitude  will  continue  to 

decrease  in  absolute  value  as  the  distance  of  the  O.P.  moves 

radially  from  the  antenna. 

The  magnitude  of  E7  and  E   are  approximately  equal 

at  the  top  of  the  monopole,  but  E   maintains  a  value  within 

an  order  of  magnitude  of  its  peak  value  all  along  the  antenna 

surface  (Figure  7).   As  a  result  of  boundary  conditions  on  E 

for  the  ground  Diane  (similar  to  those  for  E^  on  the  antenna 

z 

surface),  the  magnitude  of  E.  quickly  approaches  zero  at  the 

base.   The  shape  of  this  curve  remains  pretty  much  the  same 

as  observations  are  made  at  greater  distances.   The  magnitude 

of  the  curve  decreases  slowly  and  approaches  a  uniform  value 

along  the  antenna. 

E    .  is  a  composite  of  E   and  E  .   Near  the  antenna, 
peak         -  v  z 

Dlots  of  E    .  (Figure  8)  stronglv  resemble  plots  of  E  .   This 
peak     3  to  J  r 

implies  that  E   is  the  dominant  component  in  the  near-field 
r 

region.   However,  as  the  O.P.'s  are  moved  further  outward, 

the  E    ,  /E   resemblence  fades  and  there  is  a  pcint  around 
peak   r 

X  /  2tt  (see  Figure  9)  where  the  dominance  begins  to  shift  to 

E  .   When  E   becomes  completely  dominant,  the  value  of  E    . 
z  z  f  j  •>  peak 

is  approaching  a  uniform  value  along  the  antenna.   This 
results  from  the  electromagnetic  waves,  generated  by  the 
point  source  at  the  top  of  the  antenna,  approaching  a  plane 
wave  front  at  large  distances .   Also  noted  in  Figure  9  is  the 
sharp  decline  of  E^  above  the  height  of  the   antenna.   This 
indicates  very  little  radiation  in  the  axial  direction,  the 
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greatest  percentage  is  directed  radially  outward.   A  more 

complete  set  of  E(Z)  Plots  are  found  in  Figure  A-l  through 

A-18  in  Appendix  A. 

2.   £(R)  Plots 

Figure  10  shows  the  very  large  value  of  E   found  at 

the  top  of  the  antenna  decreases  rapidly  in  magnitude  as  the 

O.P.  moves  radially  away  from  the  monopole.   At  a  radial 

distance  of  approximately  \/.2irthe  rapid  decline  in  E^  sub- 

sides  and  approximates  the  attenuation  with  distance 

observed  at  other  points  along  the  antenna  (see  Figure  11). 

Since  E^  =  0  at  the  surface,  as  shown  in  Figure  11,  the 

magnitude  of  the  field  rises  rapidly  with  radial  distance  to 

a  peak  at  about  2  0  meters  from  the  antenna,  from  which  the 

magnitude  slowly  declines. 

The  rate  of  decline  of  E^(R)  appears  to  be  uniform 

for  all  positions  along  the  antenna,  the  only  difference 

being  the  magnitude  from  which  the  decline  begins.  (Figure  12) 

Close  to  the  antenna,  the  E    .  curve  follows  the  E 

peax  r 

curve.   At  some  unknown  point,  E  (R)  has  diminished  in 

z 

magnitude  less  than  E  (R);  hence,  the  E    ,  curve  switches  to 
&  r  peak 

follow  E  (R)  (Fisure  13).   The  point  at  which  this  switch 
z 

occurs  falls  at  a  greater  radial  distance  for  O.P.'s  at 
greater  heights.   This  is  the  result  of  a  larger  starting 

magnitude  (at  greater  heights)  from  which  E   begins  its 
uniform  rate  of  decline.   For  lamer  starting  magntudes ,  a 
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Plot  Ez  v«  r  for  a  modal  of  15  t«gm*ntt. 
Plot  !«  of  1/2  tho  monopoly  holght  oiong  tho  r  axis. 
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larger  radial  distance  is  allowed  for  it  to  decline  before 

the  magnitude  of  E   exceeds  E  .   The  ooint  at  which  This 

z  r 

occurs  is  so  sharplv  defined,  that  the  curve  of  E    ,  has  an 

-  '         '  peak 

"apparent  discontinuity"  in  it. 

At  twice  the  monopole  height,  the  curve  of  E 

to   '  peak 

(Figure  14)  has  an  inverted  cusp.   This  results  from  the  same 

shift  in  E^/E   dominance  discussed  above.   Directly  above  the 

antenna,  E   =0;  hence,  E   is  larger  than  E  .   But  as  the  O.P. 
r  z  r 

moves  radially,  E   declines  while  E   rises  to  a  oeak.   This 

forms  the  inverted  cusp;  then,  E^  also  begins  a  decline,  but 

not  as  ratid  as  E  .   More  plots  of  E(R)  are  found  in  Figures 

z 

A-19  through  A-33  in  Appendix  A. 
3.   H(Z)  Plots 

At  the  surface  of  the  monopole,  the  magnitude  of  EL 

is  virtually  constant;  there  is  a  slight  decrease  near  the 
top  (Figure  15).   Above  the  monopole,  there  is  a  sharp  drop 
(about  3  0  dB ' s )  in  intensity  of  H  ;  the  field  is  virtually 
non-existent  there.   Similar  plots  of  H.  for  greater  distances 
from  the  antenna  follow  the  same  type  of  curve,  except,  the 
magnitude  of  the  field  parallel  to  the  antenna  decreases  with 
increasing  distances  from  the  antenna,  and  the  sharp  decline 
in  intensity  above  the  height  of  the  antenna  becomes  less 
abrupt.   At  a  distance  of  one  wavelength,  the  difference 
between  the  magnitude  of  K.  oarallel  to  the  antenna  and  above 
its  height  is  about  2  0  dB.   More  plots  of  H(Z)  are  found  in 
Figure  A- 3  4.  through  A-39  in  Appendix  A. 
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4.   H(R)  Plots 

Again,  at  the  surface  of  the  monopole  the  magnitude 

of  H  is  virtually  constant.   The  magnitude  of  H  decreases 

rapidly  with  increasing  radial  distance,  and  uniformly  for 

various  antenna  heights.   Figure  16  is  an  example  curve  which 

delineates  the  behavior  of  H. all  along  the  axis.   As  indicated 

(p 

above,  for  0.?.fs  directly  above  the  antenna,  H  is  virtually 

9 

non-existent.   As  the  0.?.  moves  radially  outward,  at  twice 
the  height  of  the  antenna,  H  rises  quickly  to  a  peak,  then 
slowly  declines  in  magnitude.   More  plots  of  H(R)  are 
found  in  Figure  A-M-0  through  A-4M-  in  Appendix  A. 
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V.   A  VALIDATION  OF  NEC  NEAR-FIELD  CALCULATIONS 

NEC  is  still  evolving.   On-going  efforts  toward  refining 
and  improving  its  calculations  are  constantly  occurring.   In 
its  present  form,  it  has  been  extensively  used,  and  docu- 
mented to  provide  accurate  results  in  the  far-field;  its 
validity  in  close  to  rhe  antenna  is  not  so  well  documented. 
To  show  that  NEC  provides  reasonable  results  in  near-field 
regions,  three  tests  have  been  run.   They  are  described  in 
the  following  paragraphs . 

A.   COMPARISON  OF  NEC  WITH  CLASSICAL  ANALYSIS  RESULTS 

Mills  [Ref.  8],  using  equations  similar  to  those  found 
in  Ref.  4  (pp.  333  to  338),  calculated  theoretical  values 
for  E   as  a  function  of  radial  distance  from  the  antenna, 
for  short  antennas  of  various  lengths.   Two  of  his  plots  are 
of  interest  to  this  study;  they  are  for  dipoles  of  half- 
length  (h)  equal  to  .IX  and  ,2\.      The  results  of  his 
calculations  are  plotted  as  broken  lines  in  Figure  17,  the 
■plot  of  NEC's  calculations  for  the  same  antennas  are  shown 
as  solid  lines.   Remembering,  the  assumption  of  sinusoidal 
current  is  accurate  only  for  the  quarter-wavelength  monopole", 


"As  the  antenna  becomes  shorter,  the  sinusoidal  peak  occurs 
further  back  on  the  feed  line  and  rhe  distribution  on  the 
antenna  becomes  more  linear. 
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Figure    17 
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it  is  interesting  to  note  NEC's  calculations  are  more  closely 

matched  as  the  antenna  length  approaches  the  resonant 

condition,  especially  for  radial  distances  greater  than 

a/2tt  (=  300/2  =  48  meters  for  this  study).   Balzano  [Ref.  7] 

points  out,  that  inside  this  distance,  the  classical  approach 

breaks  down,  "...because  nowhere  in  the  near  field  (p<X/2ir) 

does  the  dipole  look  like  an  elementary  source  (infinitesimal 

dimension)  from  the  observation  point".   Balzano  also  points  - 

out  that  the  classical  approach  has  a  constantly  increasing 

E   field  ud  to  the  diDole  surface  which  is  incongruent  with 
z  ° 

energy  density  boundary  conditions  at  the  antenna  edge. 
Therefore,  the  leveling  off  of  the  magnitude  of  E  ,  as 
predicted  by  NEC,  is  the  more  reasonable  plot. 

B.   COMPARISON  OF  NEC  WITH  EMPIRICAL  ANALYSIS  RESULTS 

Balzano  in  Ref.  7  (Part  II)  performed  some  measurements 
on  a  dipole  constructed  from  bronze  rods  of  .0015  meters 
diameter  (=  3.75  X  10~  X)  that  was  .18  3  meters  in  length 
(=  .  M-712X),  including  a  .0018  meter  gap  at  The  center  from 
which  it  was  fed. 

A  model  of  this  antenna  was  constructed  using  NEC.   The 
configuration  chosen  was  a  monopole  of  19  segments,  each 
segment  being  .0  049  meters  long.   The  monopole  was  mounted 
above  a  perfectly  conducting  ground  plane  in  a  manner 

similar  to  that  used  for  this  studv.   Plots  of  E   and  E   were 

z       r 

run  as  a  function  of  Z  for  radial  distances  of  .003,  .006, 
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.01,  .02,  .03,  .04,  .05,  and  .06  meters  from  the  dipole  axis. 
The  results  of  these  plots  are  shown  in  Appendix  B.   With 
the  exception  of  a  scaling  factor,  the  information  contained 
in  these  plots  is  the  same  as  that  discussed  in  Section  IV. C. 
for  the  study.   To  further  compare  the  results  of  NEC  with 
the  measured  values,  the  measured  peaks  of  the  square  of  the 
electric  field  intensity,  for  both  axial  and  radial  polariza- 
tions, at  the  radial  distances  indicated  above,  were  extracted 
from  plots  in  Ref.  7.   These  values  were  plotted  in  Figures 
18  and  19  along  with  the  square  of  the  teak  of  the  fields 
extracted  from  the  plots  of  the  model  run  by  NEC.   As  can  be 
seen,  very  close  correlation  between  the  model  and  the 
measurements  were  obtained.   The  largest  discrepancy  was  at 
.00  3  meters  distance  where  the  calculated  value  was  inordinately 
higher  than  the  measured  value,  relative  to  other  calculated/ 
measured  pairs.   Balzano's  calculated  value  was  also  quite  a 
bit  higher,  which  would  indicate  the  measured  value  was  a  bit 
low.   Credence  for  this  argument  is  supported  by  Balzano's 
statement,  "...at  closer  distances,  capacitive  effects 
became  predominant" . 

C.   COMPARISON  OF  NEC  WITH  MAXWELL'S  EQUATIONS 

Chang   [Ref.  6]  points  out  that  "...on  the  surface  of  the 
monopole  the  field  comDonent  HA  should  be  proportional  to  the 
current  distribution  I(Z)...".   The  NEC  calculated  value  of 
HA  was  comDared  to  the  current  flowing;  on  the  model.   The 
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results  are  shown  in  Table  I. 

TABLE  I 

COMPARISON  OF  MEG  COMPUTED  H-FIELD  AND  CURRENT  AT  THE 

MO NO POLE  SURFACE 


H 


2.49 

2.573 

7.50 

2.577 

12.51 

2.552 

17.52 

2.497 

22.50 

2.413 

27.51 

2.300 

3  2.52 

2.161 

37.53 

1.9  97 

42.  54 

1.310 

47.52 

1.601 

52.53 

1.369 

57. 5^ 

1.118 

62.55 

.3478 

67.56 

.5562 

72.54 

.2247 

27rrH 


I(Z) 


4.350 

4.849 

4.858 

4.853 

4.810 

4.311 

4.707 

4.707 

4.548 

4.548 

4.335 

4.337 

4.073 

4.076 

3.764 

3.767 

3.412 

3.414 

3.018 

3.019 

2.581 

2.584 

2.107 

2.112 

1.598 

1.604 

1.048 

1.056 

.4235 

.430 

Quite  good  agreement  is  maintained  all  along  the  antenna 
The  largest  variation  occurs  at  the  end  of  the  antenna  and 
that  variation  is  still  less  than  1%. 

Chang  uses  the  difference  between: 


o  a 


Ep(p  ,z)  =  j  -:-n4iP,z 


H  .  (  P  ,  z  ) 


(63) 


and 


1  3_ 
o  o  p 


^  p  i^p  vp  ,  Z  ) 


-  iu  v.  p  ,  Z  J 


ro] 


:^a 


(64) 
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to  approximate  E  .   This  relationship  is 


yp'V  -  i  ?g  (^HW  -  H*  (V!  (55) 


Its  application  to  the  data  from  the  study  at  sample  points 
yields  Table  II. 

TABLE  II 

COMPARISON  OF  E   ESTIMATED  AND  E    CALCULATED  3Y  NEC 

P  P 


Z2 


2.5 
7.5 
12.5 
17.5 
2  2.5 
27.5 
32.5 
37.5 
1+2.5 
1+7.5 

5  2.5 
57.5 

6  2.5 
67.5 
72.  5 


$   3    *   1 

E  Est. 

0 

E   Calc. 

0 

0. 

0 

10  5 

0. 

0 

177 

.0  01 

90 

195 

.003 

270 

274 

.004 

360 

366 

.005 

450 

4  58 

.006 

540 

548 

.007 

630 

634 

.007 

630 

716 

.009 

810 

794 

.010 

900 

8  65 

.010 

300 

935 

.0112 

1008 

1006 

.0121 

1089 

1090 

.0156 

1404 

1402 

The  correlation  between  the  estimated  value  and  the  calculated 
value  is  quite  good  above  approximately  15  meters  in  height 
on  the  monopole.   The  wide  variation  below  that  point  can  be 
attributed  to  insufficient  significant  figures  in  the  values 
of  H.  used  to  estimate  the  difference.   Also,  NEC  was 
computing  the  field  emitted  from  the  base  of  the  monooole, 
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which  is  net  considered  in  the  estimate  of  Maxwell's 

equations.   Considering  these  two  factors,  the  correlation 

is  surprisingly  good. 

E   is  expected,  by  Maxwell's  equations,  to  be  zero  on  the 

surface  of  the  monopole.   In  NEC,  this  condition  is  only 

enforced  at  the  match  points.   The  best  estimate  of  how  well 

NEC  aDDroximates  E  =0  can  be  seen  from  the  plot  of  E  (Z)  at 

z  -         z 

the  surface  of  the  monopole,  Figure  20.   The  magnitude  of 

E   that  is  radiated  can  be  found  at  the  top  of  the  monoDole , 
z 

as  expected.   This  value,  from  Figure  20,  is  nominally  1000 
v/m:  at  other  points  along  the  antenna,  the  nominal  value  can 
be  approximated  as  1  v/m.   This  drop  of  three  orders  of 
magnitude,  or  60  dB,  is  a  reasonable  approximation  of  zero 
by  virtually  any  engineering  standard. 
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Plot  Ex  vl  Z  for  a  modal  of  15  togmvntt. 
Wet  It  of  on*  radius,  tho  Surfac*  of  tho  Monopolo. 
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Figure  2  0 
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VI.   CONCLUSIONS  AND  RECOMMENDATIONS 

From  the  results  of  Section  V,  it  can  be  seen  that  NEC  is 
an  effective  analysis  tool  for  near-field  calculations  and 
can  be  used  and  trusted  to  give  meaningful  results.   The 
tool  cannot  be  applied  thoughtlessly;  care  must  be  exercised 
to  ensure  the  model  accurately  represents  the  actual  physical 
structure.   Zoning  considerations  must  be  factored  into  the 
procedure  to  ensure  that  sufficiently  small  segments  are  used 
to  get  the  desired  resolution,  but  not  so  many  as  to  be 
wasteful  of  resources.   Electrical  parameters  must  be  con- 
sidered to  ensure  they  represent  the  operational  environment 
into  which  the  structure  will  be  immersea.   When  these  items 
are  adequately  considered,  NEC  will  provide  the  desired 
results  in  a  timely  and  cost  effective  manner. 

Areas  of  further  study  could  include  expanded  investi- 
gations into  the  maximum  segment  length  for  near-field 
calculations,  investigation  of  the  characteristics  of  the 
field  emitted  from  the  base  Insulator  region,  and  effects 
on  the  near-field  of  using  multiple  parallel  wires  with  cross 
connections  in  the  manner  in  which  broadcast  towers  are 
actually  constructed. 
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Plot  Ex  v«  Z  for  a  modal  of  15  t«9mont». 
Plot  l«  at  on*  rodlut.  thm  Surface  of  tho  Monopoio. 
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Figure  A-l 

Figures  A-l  through  A-6  are  plots  of  E^  as  a  function  of  z 
for  various  distances  from  the  Z^axis)V 
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Plot  Ex  v»  Z  for  a  mod*)  of  15  t*gm*ntf. 
Plot  ft  al  two  radii.  .3  M*t*ri  abov*  th*  Surfaci 
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Figure   A- 2 
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Plot  Ez  vt  Z  for  o  mod«i  at  15  »«gm«nti. 
Plot  It  at  flvo  radii.  1.2  Motvrt  abov*  th*  Surfaco. 
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Plot  Ex  vt  Z  for  a  mod»l  of  15  togmontt. 
Plot  it  of  t«n  radii.  2.7  Motors  abovo  tho  Surfoco. 
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Figure    A-M- 
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?k»t  El  v*  2  for  o  ato&i  of  15  t«^a«ftH, 
Pto«  h  o»  tOO  rwin.  29 J  Motor*  obovo  ?*•  Swrfooo. 
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Figure    A-5 
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Wot  &viZ^r«  n*d*l  «tf  15  ««ge*«n*«. 
Plot  H  at  WOO  rvdn.  299.7  M«t«n  obov«  rf»  3ur««c» 
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Figure  A-6 
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Plof  Er  y«  Z  for  o  modol  of  15  togmonM. 
Plot  I*  at  on*  rad!u«.  tho  Surfaco  of  tho  Monopoly. 
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Figure   A-7 

Figures  A-7  through  A-12  are  plots  of  E_  as  a  function  of  z 
for  various  distances  from  the  Z-axis. 


94 


Plot  Er  vi  Z  for  a  modol  of  15  sogmvnt*. 
Plot  It  al  two  radii.  .3  M«t«n  abovo  th«  Surfaco. 
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Figure   A-i 
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Plot  Er  v*  Z  for  a  mod*l  of  15  sogmonts. 
Plot  It  at  flvo  radii.  1.2  Motors  abovo  tho  Surfaci 
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Figure   A-9 
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Plot  £r  vi  Z  for  a  modol  of  15  togmontt. 
Wot  It  at  ton  radii.  2.7  Motor*  abovo  (ho  Surfaeo. 
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Figure   A-10 
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Figure   A- 11 
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Plot  .»  «t  !OC0  reoU  299.7  Motor*  «Jx>v»  tho  Wfac* 
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Plot  Epook  vt  Z  for  a  medoi  of  IS  i»gni#nti. 
Plot  It  at  on*  radlui.  th«  Surfac*  of  fh«  Monopolo. 
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Figure   A-13 


Figures  A-13  through  A-13  are  plots  of 


me  ax 
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of  z  for  various  distances  from  the  Z-axis. 
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Plot  Epook  vt  Z  for  a  modal  of  15  tvgmants. 
Plot  It  at  two  radii.  .3  M«tors  abovo  th«  Surfac*. 
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Plot  Ep«ck  v«  Z  for  a  mod«l  of  15  i*gm»nt». 
Plot  It  at  flv«  radii.  1.2  M*f»r«  abov*  th«  Surfac*. 
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Figure   A- 1 5 
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Plot  Ep«ok  vt  Z  for  a  mod»i  of  15  t«gm«nf t. 
Plot  It  at  f«n  radii.  2.7  M«f«r«  abov*  th*  Surface. 
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Plot  Epodk  w  Z  for  «  — dW  of  ?5  sogaontt. 
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Plot  Ez  v«  r  for  a  model  of  15  togmontt. 
Plot  Is  at  2/lOOth  tho  monopoio  holght  along  tho  r  axlt. 
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Figure  A-19 


Figures  A-19  through  A-23  are  plots  of  E^    as  a  function  of 
for  various  heights  above  the  ground  plane. 
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Plot  Ei  vi  r  for  a  modal  pf  15  sograonti. 
Plot  Is  at  1/4  tho  monopolo  holght  along  tho  r  axlt. 
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Figure   A-20 
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Plot  Ex  vt  r  for  a  mod«l  of  15  «*gm«nti. 
Plot  It  of  1/2  th«  monopolo  h«tght  aiong  th«  r  a*l». 
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Figure   A-21 
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Plot  El  vs  r  for  a  modal  of  15  ••gmonts. 
Plot  It  of  tho  holgbt  of  tho  monopolo  along  tho  r  axis. 
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Plot  Ex  vir  for  a  modal  of  15  togmont*. 
Plot  1»  of  fwleo  tho  monopolo  holght  along  tho  r  aift. 
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Figure   A-2  3 
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Plot  Er  v«  r  for  a  modol  of  15  t ogmonH. 
Plot  (•  of  2/lOOfh  tho  monopelo  h«lght  along  tho  r  axit. 
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Figure   A- 2 4 

Figures  A- 24  through  A-28  are  plots  of  E^  as  a  function  of  R 
for  various  heights  above  the  ground  plane. 
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Plot  Er  vi  r  for  a  mod«l  of  15  s»gm«ntt. 
Plot  It  of  1/4  tho  monopol*  h«ight  along  tho  r  axit. 
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'MO- 


Plot  Er  vi  r  for  a  modot  of  15  *ogm«nf«. 
Plot  1»  at  1/2  fho  monopoio  holght  along  th«  r  axlt. 
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Figure   A-25 
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Ptot  Er  v«  r  for  a  modal  of  15  t«gm«nt«. 
Plot  I»  at  tho  height  of  tho  monopoU  along  tho  r  axlt. 
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Plot  Er  v»  r  for  a  modal  of  15  t«gmont«. 
Plot  It  at  twlco  tho  monopoto  holght  along  tho  r  axis. 
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JO«H 


Plot  Epook  y%  r  for  a  mod«l  of  15  i»gm«nti. 
Plot  It  at  2/100th  tho  monopoU  height  along  tho  r  axit. 
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Figure  A-29  through  A-33  are  plots  of  E    . 
for  various  heights  above  the  ground  plane. 


as  a  function  of  R 
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Wot  Ep«ok  vt  r  for  a  modal  of  15  t«gm«nt«. 
Plot  la  at  \/A  tho  monopolo  holght  along  tho  r  oiii. 
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Plot  Epook  v«  r  for  a  modoi  of  IS  i«gmontt. 
Plot  ll  ot  1/2  fho  monopoio  holght  along  th«  r  axlt. 
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1000- 


Plot  Epook  vi  r  for  a  modal  of  15  togmont*. 
Plot  It  at  tho  holght  of  tho  monopoio  along  tho  r  axis. 
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IOt 


Plot  £p«ak  vt  r  for  o  modal  of  15  togmontf. 
Plot  It  al  twlco  tho  monopol*  holght  along  tho  r  axis. 
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Plot  HpM  vt  2  to  Twlco  onfonno  K«lght.  13  t«gm«nH. 
Plot  ft  of  ono  rodhia.  tho  Surfoc*  of  th»  Monopol*. 
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Figures  A-34  through  A-39  are  plots  of  H   as  a  function  of  Z 
for  various  distances  from  the  Z-axis . 
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Piof  Hphl  vi  Z  to  twice  ontonno  hol^ht.  15  togmont*. 
Plot  I*  at  two  rodU  .3  M«t*rt  obevo  th*  Surfoco. 
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Plot  Hphi  tilto  twlc*  entonna  height.  13  j«g*«nrt. 
Wot  It  at  flv«  radii.  L2  M«t»rt  obov*  thm  Surfc 
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Plot  HpM  v«  Z  to  twlc*  orrf «noa  tolght.  13  I •gmontt. 
Rot  li  at  ton  rodlf,  2 J  M»t«rt  cbov*  th*  Surfoc*. 
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Wot  Hpfcl  v%  2  (©  twte*  oottnnq  height.  IS  isqm+ntt. 
?!ot  li  at  WO  rodn.  29.7  M*4ec»  ■**»**  Hm  Surfoc*. 
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Ptot  HpH  r%2  f  twice  «rrt»rwi«i  height.  t3  9«g«e<%.«. 
Ptot  is  at  KXX)  radii.  299.7  Meter*  above  rite  Surface 
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Plot  Hphl  vt  r  for  a  mod*)  of  15  sogmontt. 
Plot  It  ot  2/100th  th«  monopoly  height  along  tho  r  axil. 
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Figure   A-^0 

Figures  A-40  through  A-^u  are  plots  of  H   as  a  function  of  R 
for  various  heights  above  the  ground  plane. 
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Plot  Hphl  vt  r  for  a  flto<4*l  of  13  togmonts. 
Pfot  I*  at  1/4  tho  mooopoJo  Hoight  along  tho  r  ail*. 
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Plot  Hphl  v«  r  for  o  aoool  of  15  togmontt. 
Plot  H  ot  1/2  tho  aonopoJo  height  along  tho  r  axis. 
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Wot  Hphl  v«  r  for  a  mod*)  of  15  togmonH. 
Plot  It  at  tho  holght  of  tha  monopoio  along,  tho  r  axle 
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Plot  Hphl  virfora  aodol  of  15  aognonts. 
Plot  It  at  twfco  tho  nonopoio  hoJght  along  tho  r  cod*. 
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APPENDIX    B 


PWt  flvit  for  o  — rftl  of  Y9  togaoott. 
h  .3  en  frwa  ■owpoU  odi  d»n§  riv»  s  oaf*. 
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Figure    3-1 

Figures    B-l    through    B-8    are    plots    of   ! 
for    various    distances    from   the    Z-axis 


as    a    function    of   Z 
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Plot  Eini  for  «  — 4*i  ot  19  t*9«**ti. 
Pfof  It  ^em  fro*  rh#  «»mpo)i  writ  afong  th«  i  aattt. 
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Plot  Ex  *•  x  for  a  «oo*J  of  19  ••gmantt. 
Plot  '*  t  em  freai  Hm  »onopoi«  asli  aiong  th«  s  cart*. 
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Plot  Eiyii  for  a  aocM  of  19  i«g*«ntt. 
Plot  it  2  am  from  tho  wonoaoU  axlt  along  tKo  z  axil. 
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Hot  Ex  v«  t  for  o  aodo!  of  19  sogaoott. 
Plot  ii3«a  froio  tk*  »onopol»  ost*  oiong  th»  s  oxlt. 
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Plot  n»ti  for  a  aod*4  el  19  i«g*«nt«. 
Plot  It  4  cm  fr»w  tK*  wonopoi*  «rt»  along  *ho  x  «xk. 
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Plot  ft  5  cm  1rom  thm  aonopofe  writ  qsong  tK*  i  corit. 
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Plot  brii  for  a  m*4m\  of  19  ••gaiontf. 
Plot  to  6  cm  from  tfc*  ■onopoU  «xJ«  along  tho  z  aads. 
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Plot  £r»u  for  a  m*6*l  of  t°  t+qmfttu 
It  .3  ob  fro*  rt*#  — wopolo  oadt  oiong  tho  x  osit. 
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Figure  B-9 
Figures  B-9  through  3-16  are  plots  of  E   as  a  function  of 


for  various  distances  from  the  Z-axis. 
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Plot  Er  vt  x  for  o  mo6mi  of  19  log— wto. 
Plot  It  .6am  trom  tfco  •oooooio  oils  oioog  tho  x  «xl«. 


0.00 


0.03 


1  T  ■  ~ 

0.04 


— I 

0.04 


0.00 


0.10 


o.n 


DHtonco  Z  (a) 


Figure    3-10 


141 


Plat  Er  v*  i  for  a  w—frl  of  19  fynf. 
Plat  ft  I  an  fro*  tho  ■aaofaia  «***  along  tna  i  wdu 
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Figure    3-11 
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Plot  frru  for  a  mmdaH  of  t9  foeMontt. 
Plot  it  2  an  horn  tfco  wonopoio  «n»  oiong  tfco  z  end*. 
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Hot  Er  y*  t  for  o  aodd  of  19  tagaontt. 
Plot  It  3  ca  fro*  tha  aowopofc  «jri»  otooo;  tha  z  writ. 
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Plot  Er  ¥%  i  for  a  aocfei  of  19  ••gin«nft. 
Plot  H  4  cm  fr*«i  th*  »«nopo*«  axi»  along  tho  x  axle 
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Plot  Er  v«  z  for  a  — dot  of  1°  togoiontft. 
Plot  Is  5  en  frow  tko  ojoooooJo  mdt  along  tho  z  azli. 
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Plot  Er  v«  s  for  o  iocUI  of  19  »«g«»n»». 
Plot  to  6  «■  froMi  tfto  MOAopoJo  osit  a*ong  tho  z  axH. 
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